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ABSfRAC i 



‘i study is made of compressible perfect fluid motion in 
turbomachines having infinitely many blades and a general theory 
is developed. An underlying concept of the theory is that force 
fields which represent the action of infinitely many blades belong 
to a special class described as pseudo-conservative and can be 
expressed as the product of a scalar function and the gradient of 
a potential. The scalar function is simply the rate at which energy 
is imported to the fluid by the blades, and the potential is simply 
the family of the equations for the blade surfaces. The introduc- 
tion of these two functions to express the force field casts an en- 
tirely new light on problems of mixed-flow turbomachines having 
infinitely many blades of arbitrary shape. 

In the formulation of the problem the non-linear action of 
rotationality and compressibility is regarded as a force tending 
to displace the streamsurfaces from their irrotational , incom- 
pressible position. It is shown that the character of the problem 
is determined by a governing velocity: the velocity relative to the 
blades where blades are present, or the meridional velocity, 
where blades are not presenr, , Where the governing velocity is 
subsonic the problem is essentially elliptic, where supersonic, 



hyperbolic . 



The theory and the examples lead to conclusions which 
are believed to explain in part the unexpected efficiencies ob- 
served for compressors having transonic governing velocities. 
These conclusions, which indicate that transonic compressors 
could perhaps be profitably developed, are as follows: The de- 

flection of the streamsurfaces induced by a given strength of 
vorticity at a certain point in the flow has one sense when the 
governing velocity at the point is subsonic , the opposite sense 
when it is supersonic, and becomes aero as it becomes sonic. 
The deflection of the streamsurfaces brought about by a given 
distribution of vorticity in a region is less when the governing 
velocity in the region is transonic than when it is entirely sub- 
sonic or entirely supersonic. 

Examples of incompressible flow through a mixed flow 
compressor with prescribed blades, and subsonic and transonic 
flow through actuator disks , were solved by the method of 
finite differences, applying simultaneously the relaxation 
technique and an iteration process . 
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i. iNTRODUC f . ,-N 



In a real turbomachine a viscous , compressible fluid flows 
through an axially symmetrical channel. In a region of this channe 1 
a system ox blades, either stationary or rotating about the -xis of 
symmetry, acts on the fluid. The field of the xlo'\ , being bounded 
oy the surfaces of the blades as well as by the channel boundaries, 
is not circumferentially uniform. The total energy of the fluid may 
vary from point to point in the field, the flow is genera -ly rotational, 
and if the blade system rotates, the flew is unsteady. I’ he real 
problem is thus a very formidable one, and cannot be approached 
exactly by methods known today. If viscous effects are neglected 
the problem is greatly simplified, because of this, and because 
viscous effects actually are very slight except near the boundaries, 
the assumption of zero viscosity has always been made . 

/■ brief but excellent review of the earlier important investi- 
gations of the flow in turbomachines io given by ivarble (Reference 1) . 
Earlier investigations were concerned primarily with the flow through 
a typical annulus of smalt radial extent and hence treated the flow as 
essentially two-dimensional , e. g. Traupel (Reference 2). Interfer- 
ence between the flow in neighboring annular regions was assumed 
negligible, a condition which is fulfilled only if the centrifugal iorces 
and the radial pressure gradient forces are in balance, ;ne.t is, if 
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the flow is that of a vortex centered on the axis of rotation. This 



is a rather severe restriction, too strong to be imposed in most ap- 
plications . 

The first detailed analysis of the three-dimensional incom- 
pressible flow in turbomachines was given by !v eyer (Reference 3). 
An exact solution was- obtained for machines in which the flow is ir- 
rotational upstream and downstream of the rotor, although the flow 
within the rotor may be rotational, h eyer achieved circumferential 
uniformity by assuming an infinite number of very thin blades, and 
then modified this result by a Fourier analysis to obtain the solution 
for a finite number of blades . This was an important contribution 
but the extension to arbitrary blade shape, with rotational flow down- 
stream of the rotor, has not been made. , 

Marble (References 1 and 4) linearized the rotational incom- 
pressible problem with an infinite number of blades by the assump- 
tion that the vorticity is transported along the streamlines of an ir- 
rotational flow within the same boundaries . Ke provided examples 
of axial flow and conical flow, and stated that the 3imple linearized 
solution was sufficiently accurate if the vorticity effects were not 
large. A second order linearization was given to handle problems 
in which the vorticity effects are large. Using the simple linear- 
izations Marble considered three interesting problems: the mutual 



interference of neighboring blade rows in a multistage axial turbo- 
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blades is completely arbitrary. The two assumptions -.'.hic’ri .naxe 
the analysis possible are: the fluid considered is non-viscous , and 

the number of blades is infinite. The viscosity of the fluid is unim- 
portant except within the boundary layers along all the bounding sur- 
faces , and these boundary layers are thin if pressure gradients are 
favorable. F urtliermorc, an elementary consideration of boundary 
layers often presupposes a knowledge of the velocity distributions 
such as obtained by this analysis. The assumption that the number 
of blades is infinite makes the flow field circumferentially uniform. 
This uniformity makes the solution easier in that the flow depends 
only on two coordinates, but at the same time it introduces arti- 
ficial complications . Generally there will be a discontinuity in 
the flow as it enters this blade region if the number of blades is 
infinite. Furthermore the force o.‘ the blades on the fluid is not 
applied on individual blade surfaces, but is distributed throughout 
the blade region and actually acts as a body force field distributed 
throughout the region. The circumferentially uniform flow field 
may be considered as the limiting case of a flow acted upon by a 
very large number of closely spaced blades of negligible thickness. 

Two distinct problems occur in practise. The design prob- 
lem arises when it is desired to design a machine for a particular 
purpose . i’he analysis problem arises when it is desired to 
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invcstigate a given machine under given operating conditions. These 
problems were classified by Marble {Reference 4) as the indirect and 
the direct problem, respectively, in analogy with the three-dimen- 
sional wing theory. This classification is convenient and very sig- 
nificant. Here any problem in which the blade shape is prescribed 
will be classified as a direct problem, any other problem as an in- 
direct problem. The full significance of this distinction becomes 
clear when the circumstances under which the idealized problem is 
actually comparable to a real problem are understood. 

In the real machine, which has a finite number of blades, 
each blade transmits a force to the fluid by maintaining a discontin- 
uity in pressure across its two surfaces. If no viscous force is 
present this blade force must be normal to the blade surface. In 
idealizing the problem these concentrated blade forces are essen- 
tially replaced by a body force field. Clearly this body force field 
must be normal to the family of blade surfaces throughout the region 
of the blades if the two problems arc to be comparable. It is im- 
portant to note that this imposes a purely geometrical restriction 
on the body force field. It will be shown that the necessary and suf- 
ficient condition that it be possible to construct a family of surfaces 
which are everywhere perpendicular to a given force field is that the 
force field is everywhere perpendicular to its own curl. If the blade 
shape is prescribed, as in the direct problem, this offers no difficulty. 
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However in the indirect problem where the blade shape is one of 
the things to be determined, the prescribed quantities, which for 
example, may be components of the force field or the energy dis- 
tribution, must be prescribed in such a way that the existence of 
the family of blade surfaces is assured beforehand. This necessity* 
was noted by i.!eyer (Reference 3) in connection with the prescrip- 
tion of vorticity in the indirect problem, but was overlooked by 
Reissner (Reference 6), Gravalos (Reference 7), and several others. 
When the streamwise extent of the working region is much less than 
its transverse extent, that is, if the blade aspect ratio is large, as 
in most axial flow problems , the idealised flow will be comparable 
to the real flow even though the blade existence has not been assured. 
Thus garble (Reference 4) was able to obtain good approximate 
solutions of the indirect axial flow problem by prescribing blade 
force components. However in mixed-flow problems the aspect 
ratio is not large, and in the indirect mixed-flow problem the exist- 
ence of the blades must be assured beforehand. 

In the real problem one boundary condition is that the fluid 
does not flow through the blade surfaces . This cannot be applied 
as a boundary condition in the idealized problem oince the blade 
surfaces, being infinite in number, are not boundaries of the flow . 

* Information lately received reveals that this condition for integra- 
bility of the blade surfaces v/as stated by lfauersfeid as early as 
1905 in a paper (Reference 10) not available at the time of writing. 



instead th* fluid in the region of l.V blade? is constrained to r -,,r - 
‘ rient 5a thc> biad * surfaces. This catmint, quite different from 
toe boundary condition of the real problem, plays an important part 
in the mathematical formulation of the ide?.Li?,ed prot'e. .. The 
geometrical relation ••’hicu expresses the constraint is equal ir 
significnnce to the ether equations •°nd must he c or s id. ’rod simul- 
taneous' y with them. In addition, ,-ith uniform inlot esnditiour and 
'- n * nc ornpre:? # il>le fluid, the differential equation for the flow lr 
'inear in the region of the blades, although it in generally non- linear 
Uov/nstre i.v of this regicr. . This fact a, v tied ’ 'oyer (:? efermcn 3) 
tc find exact solution* fur flow* in w.hic* tbs wrficity dev rstres-n 
oc the blades ir zero, although thr flrv .may rotational nud very 
complicated in the region of the blades. 

To su'..-iT ; t r :ri?* , the idea lie 2 d problem must satisfy the equa- 
tion* of motion with the b force field included, the continuity 
aquatic n, the u entropic prcaau.m-density relation, and these bnun- 
d'.ry conditions ci the roa! problem not applied on tho blade surfaces. 
il.tr two additional requirement? which must be fulfilled if the so- 
lution of the idealised problem is t.> bs co r.pvrnble tc. the real prot«- 
l«.v. are: the lody force field vas- be ncr~r.nl to th. family of blade 

surfaces; and the velocity vector field must he parallel to the in ro- 
ily tn blade surface. a. 

in the direct problem the blade vh pa ia kuo va uni -h • 
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supplementary conditions can easily be irrposcd. Por the indirect 
problem in which force components are prescribed these require- 
ments should be expressed in terms of the force components as 
follows: the body force field must be normal to the velocity vector 
field; and the force vector must be normal to its own curl. These 
assure the existence of a family of surfaces vhicli can be chosen 
to fulfill the two requirements stated above for the direct problem. 

If in the indirect problem ether quantities such as energy distribu- 
tion are prescribed, equivalent requirements must be satisfied. 

The idealised problem is formulated in terms of a stream 
function for the %'elocities in the meridional plane, and the result- 
ing differential equation for the stream function is written as a 
second order non-homogeneous partial differential equation by re- 
garding the non-linear terms, representing the effects of rotationality 
and compressibility, as forcing functions tending to displace the 
streamsurfaces . The differential equation is then replaced by a 
finite difference equation which is solved by a simultaneous appli- 
cation of the relaxation technique c£ Southwell (Reference 8) and an 
iteration process . 

The general theory is developed in Parts II through VI, and 
the difference formulation and examples are presented in Parts 
VII, VIII, IX, and X. The examples were conceived in order to 
demonstrate separately the different phases of the problem. Part 
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VIII is concerned with incompressible flow through a region of 
prescribed blades and is primarily an example of rotationality 
in a mixed-flow compressor. Part IX is offered as an example 
of the effects of compressibility when the vorticity distribution 
is approximately constant. Part X is an example of the inter- 
action of rotationality and compressibility effects when the total 
velocity is transonic but the governing velocity is subsonic. 
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ihe flow is described (Figure 1) for the most part in an ab- 
solute cylindrical coordinate system r, 6 , z by the velocity compon- 
ents u, v, v/ respectively . The absolute velocity vector is: 

[/ = i t U +■ 7 e V *- t z W , 

where i r t i 9) and i z are unit vectors in the r, 0 , and a directions 
respectively. The vorticity components for circumferentially uni- 
form flow are respectively: 



Figure 1 

The Velocity and the Vorticity Components in the Absolute Coordinate 

System 

The scalar functions u, v, w, 5 are independent of the 0 

coordinate in accordance with the requirement of circumferential uni- 
formity . 



~ • n - <Dw (■ _ / 2> rv 

~ Oz Or j J r Dr 



( 1 ) 



and the vorticity vector, _Q = v*V , is: 



n = t r / v- 4 7 + 7 Z | 
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The forca of the blades on the fluid is given by 

F ~ T r F r +7 e Fg +■ T z F z 



and has the dimensions of force/unit mass, or 



length 



time** 

At times the relative cylindrical coordinates r, z, where 
<j> is measured with respect to the rotor, will be used. If the angular 
velocity of the rotor is u) , the absolute and relative velocities are 
related by : 



do - c/ <p = U> di (2) 

The velocity and" vorticity vectors in the relative system are respec- 
tively: 

(A/ = V - 7+ cor = T r a + (v- u)rj + 7 Z w J and 

fl ’ jTL-UZ oO = T r £ + 7 j q + 7 Z - £ cC ) 



The intrinsic coordinates n, 8, s will often permit a more 
concise presentation (Figure 2). n is distance normal to the merid- 
ional streamline and s is distance along the miridional streamline. 
The velocity and vorticity vectors are respectively: 

V = iff ^ + T s q s } and 

-o - -O n + 7 B >2 + U sis j 



where 

-^n 



7 



j_ dry 
r t*s j 



- 



~ *s 9s 



^<75 

dn 



j_ dry 
r O n > 



( 3 ) 



and where K g is the curvature of the meridional streamline. 

The family of the blade surfaces is given by a scalar function 






i 
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Thc Velocity and the Vorticity Components in the Intrinsic Coordinate 

System 

of position, called the blade surface function: 

A blade surface is defined by {3 = constant. The function f(r , 2 )= (0-/3) 
defines the family of the traces of the blades in a meridional plane. 

The normal to the blade surface is 

V /3 ~ - T r F r + ~ 9 ~r ~ '* 

Since the blade force F is always normal to the instantaneous sur- 



face of the blades it can be expressed as 

F - X(r,z) V/3(r,e,z) ( 6 ) 

This defines the scalar function /^(r,s), From the second equation 
of motion it will be shown later that cOA is the local rate at which 
rey is added to the fluid by the action of the blade forces. (3 



energy 



and f are non-dimensional, 7(3 has dimensions of (length) , and 



h' a dimension of 



length'' 
time 4 ' 

; v ote that vfi , the normal to the blade surfa ct , i.-> nut i. ufeit 
vector . Instead it is chosen so that its tangential co.:k orient ij ; -.it . 
Thus the magnitude of the tangential component oi the force vector 
is simply F& - yf A . 



1 iscellanaous ymbois : 



Z 

P 



r 



if 
Z r 



if 

tu 



- static pressure 



p = density 



a = speed of sound 

, = ratio of specific heats 

t = time 

= - total energy of the fluid 

( ) j subscript denotes condition far upstream of the region 
oi the blades, or at the upstream boundary 

( ) t subscript denotes condition at the trailing edge of the 
blades 



( ) r ( ) * partial differentiation with respect to r and a re- 

spectively, unless otherwise defined 

a = angle between the blade surface and the meridional 
plane measured in a plane where s is constant 

a n = angle between the blade surface and the meridional 

plane *• leasured in a plane tangent to the surface vhoru 
r it cunst-i r t 

= angle betv/eeu the blade surface and the meridional 
plane me'-sured in a plane normal to their lim. of ir.V’r- 
section 
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= curvature of the meridiem i line 

£ 

K - curvature of the normal to the meridional street aline 

1 “ meridional velocity component 

& 

^ = ,l residual’ of the difference operator 

.^( ) ~ difference operator 

r = circulation about an annular vortex ring 

( T = angle between meridional velocity and axis 

a.. = influence coefficients 

M = Mach number 

r = outer radius 
o 

r^ = inner radius 

r = — — r, ~ - non-dimensional radius 

Hz ” 
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m. n. - fseudc-ogku^rvativ v force field 

' ND THE SLADE 3U Ax' AC J FUNCTION 

.A) The Necessary and .Sufficient Condition for the Existence of the 
arr.ily of Blade surfaces 

’’he body force field F is normal to the family of blade sur- 
faces /3 , hence F and vfi are parallel and can be related* by 
^ ~ ^( r j z ) 7'/3(r J & J z) 

where /A(r,z) is a scalar function, and vfS> is a function of only r 
and z because of circumferential uniformity, 'writing 

lirjF « 7/3 

and taking the curl of both sides there results: 

Or) V x F + * (it)* F - d? 

I* 

Then forming the scalar products with F , 

(wJ F- v * F = O 

and since A is. not zero unless |F| is zero. 

F ■' r A F = o ( 7 ) 

This is the necessary condition for the existence of the family of 
blade surfaces and is the additional relation which must be fulfilled 
in the indirect problem if the idealised problem is to be comparable 
to a real problem. 

In order to discuss the question of blade existence suppose 
that in the indirect problem force components are to be prescribed. 
First consider flow through the channel with no force presen*, fhe 

fl_o\/_is completely determined by the equations of motion, the con- 

* The applicability of Equation (6) and the subsequent derivation of 
Equation (7) were pointed out by Dr. Arthur Erd&lyi in personal 
conversation . 






*1 
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tinuity equation, a pressure-density relation, and the necessary 
boundary conditions . Now consider flow through the same channel 
with forces present. Three unknowns, the three force components, 
have been added to the problem. One equation, the requirement 
that the relative streamline be normal to the force field, has been 
added: 



At this point it appears that two components of the force can be 
prescribed arbitrarily. The third component would then be expressed 
by means of Equation (8), and the problem would be completely de- 
termined as in the channel with no force acting. The solution would 
describe a flow under the action of the prescribed forces, but would 
not necessarily be comparable to a flow acted upon by blades because 
the existence cf the blade surfaces has not been assured. If Equation 
(7) is imposed the problem is redundant, for with two force compon- 
ents prescribed there are more equations than unknowns. An example 
will give physical meaning to this difficulty. Suppose the two com- 
ponents F r and Fg are prescribed as: 



F k ® ) F-p — F& ( 

F z can be expressed in terms of these and the velocity by Equation 
(8) as: 



W * F = o 



( 8 ) 



jfk. 

F 9 



v — cd r 



(9) 




I 
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fhe three force components are then known and the problem is com- 
pletely prescribed as in the channel with no force acting. The solu- 
tion will give u, v, w as functions of r and a. Now the equation as to 
the existence of blade surfaces is raised. Equation (7) for this case is 



>\ r r d r ~~ 7 



dr 



from which 






= function of z only 



( 10 ) 



The problem is clearly over prescribed, for F a is restricted when 

Equation (10) is imposed. Furthermore, it is clear what should be 

prescribed. If P r is zero, then —~r can be a prescribed function 

rrQ 

of z only. From Equations (9) and (10), 



r 



- -r x function of z 



* fe 

This function of z that may be prescribed represents an angle be- 
tween the relative velocity and the meridional plane and actua ly is 
the blade shape parameter f z (z) defined below. Thus when the exist- 
ence of blades is assured the indirect problem with force components 
properly prescribed is equivalent to the direct problem with the 
blade shape prescribed. 

B) The Viost General Form of the Body Fo rce Field 

The reason that the two force components cannot be indepen- 
dently prescribed is that the force field, although non-conservative, 
is restricted in form. The most general force field that will satisfy 
Equation (7) is given by Equation (6). This is the most general force 
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field for which a family of blade surfaces actually exists and is 
therefore the most general force field possible in turbomachine 
problems where the number of blades is infinite . Tor this reason 
it is proposed for all such problems that the body force field be 
expressed as the product of a scalar function and the gradient of 
a potential: 

S = A v/3 = \ v [s -f (r ,?)] ( 11 ) 

This force field has a special place between conservative fields 
and non-conservative fields. conservative force field can be 
expressed in terms of one scalar function, its potential. A non- 
conservative force field requires three scalar functions for expres- 
sion, one for each component. 1 he force field of the turbomachine, 
being intermediate between these in that it requires two scalar 
functions for expression, might be called a pseudo-conservative 
force field. Actually its vectors have the same direction as a con- 
servative field, but not the sarnie magnitude. 

This fact that the force field depends on only two scalar 
functions indicates why , in theindirect problem with force compon- 
ents prescribed, the two components cannot be prescribed indepen- 
dently. It also indicates that considerable simplification might 
result from considering the functions X and f instead of the three 
force components. This is actually the case. The introduction of 
these two functions for the force vector field casts an entirely nr w 
light on turbomachine problems in which the number of blades is 
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infimte . For exarnp.e , the heretofore more difficult direct prob- 
lem becomes relatively simple, furthermore the functions A and 
f have special physical significance: c^A(r , z ) represents the local 

time rate at which energy is added to the fluid by the. action of ro- 
tating blades, and f(r,z) describes the blade shape completely. 

1' rom Equation (11) the force field is finally expressed as 

~ ^ (~ *r tr + 'e ~pr / 2 J ( 12 ) 

\vhen the force field is expressed in terms of the functions 

A and f the question of blade existence dees not arise and the direct 
and indirect problems become quite clear. / direct problem is one 
m which the blade shape f(r,z) is prescribed and X(r ,z) is an un- 
known determined by the second equation of motion after the solution 
is complete. The function A does not appear in the equation for the 
meridional flow. An indirect problem may be one in which the local 
rate of energy input, A(r,z) is prescribed, and the blade surface 
function, f(r,z), is an unknown dependent variable . 

— ) fhe ?h y a >cal leaning of the Blade Surface Function /ft ; 

It is of interest to note the geometrical significance of {3 
and f and their derivatives. The absolute value of v/3 is 

Jy/sJ - -£-[ /+.#* + Z 2jd- ( 13 ) 

where /? = r f r and Z = 

i he relation between the vector v/3 and the blade surface is shown 

in Figure 3. Using the definitions of a o , a as given in Part 

A z n 
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II it easily follows that 



R = tan c 

r 

Z = -tan a 

z 

%vith the sign convention shown 



r 




Looking Downstream onto 
a Plane z = constant 



04 ) 



re 




cooking Radially Inward onto a 
Plane Tangent to a Surface 
r = constant 

3 



Figure 

The Slade Surface Function f(r,z) 
In the same way 



R 2 + Z 2 = -/-an 2 °< n (15) 

and from Equation (13) 

r / 7 ft] - sec. r> 0 ^) 

Thus R and Z are merely the tangents of certain angles be- 
tween the blade surface and the meridional plane. 



The functions /3 and f were defined in such a way that /3 = 



constant iG the equation of a blade surface and f = constant ic the 
equation of the trace of a blade surface in the meridional plane . 



I 
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If the blades are so-called radial blades, that is, if the blade sur- 
faces are generated by lines that are normal to the axis, the function 
f is independent of r . In this case f r = 0 and f 2 is a function of a 
only . 



I 
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IV. THS HYDRODYNAMIC • L ZC U^I'ICNS 



-\) The equations in the Various Coor dinate bystems 

The hydrodynamical equations for the steady, adiabatic, 
axially symmetr ical , and circumferentially uniform flow of a non- 
viscous fluid acted upon by a pseudo-conservative* body force field 
are given in the various coordinate systems previously described. 
Vector Forms: 



The equation of motion: 

V ■ v V - - js- vp + A r/3 (i?a) 

V * -Q = P- 1/ + vp - A 7/3 (17b) 

The continuity equation: 

v • t° V - O (18) 

The isentropy condition: 

<’(—,} =0: a 2 = ££ ‘ ll-'ljf - p'" (19) 

Absolute Cylindrical (bordinates: 



The equations of motion: 



u 



u 



u 



d U t> U yZ. 


/ <2/0 


AS r 


— + w — - — = 


“ 7~ Dr 


dr d Z 


X 




d w 2) W 


/ 


A f x 




“ r> Dz 



(20a) 

(20b) 

(20c) 



The continuity equation: 

-ri?( pru ) + = ° (2i) 

•fore*. 

* a pseudo- conservative freisoo field is defined as one which can be 
expressed as the product of a scalar function and the gradient of a 
potential, lee Part III. 
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•\osolute Intrinsic (F low) Coordinates: 

The equations of motion: 

9s - ~T~ sin <T = ~T fa 

xv 2) rv v 

9s 2 s ~ ^ 

K s qrf - -j? cos <T - ~ J *2% ~ 

The continuity equation: 

Js ^ >efs + H„ t- P q s -jrs)n<r = O 
B) An Integral of the E quation of Motion 

uince the flow is steady the equation of motion can be inte- 
grated along each streamline. Forming the scalar product of V and 
Equation (17b), a scalar equation is obtained: 

V- 7 TVV- f -~p-V-vp-AVr/3 = O (24) 

From Equations (6) and (8) the condition that the relative velocity 
vector and the force vector are perpendicular is 

W - 7/3 - ( V- 7 0 u)r ) ■ v /3 = o ( 25 ) 

Consequently 

V V{3 = } and v- u>r = rP r u + (26) 

From the second equation of motion (Equation (20b)): 

^ = * 7 '- *(rv) < 27 ) 

Since the flow is isentropic the density is a function of pressure only, 
p = p(p), and the pressure term of Equation (24) is: 




(22a) 

(22b) 

(22c) 

(23) 
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Using Equation! (26), (27), and (28), Equation (24) becomes 

V ■ V [ +- y z + jj*£L _ a) rv] = O ( 29 ) 

Integrating along the meridional streamlines an energy equa- 
tion analogous to the Bernoulli equation is obtained. 

■jr V * J -~£ - cd rv - constant on each strearnsurface (30) 
Upstream and downstream of the blades rv is constant on each stream - 
surface . The total energy is E = . rhe term V s is the 

kinetic energy per unit mass and is the enthalpy. The change 

of <Jrv which occurs along a strearnsurface in the region where 
blades act is the change in the total energy of the fluid occurring on 
the same strearnsurface * 



j. he constant in Equation (30) will have the same value on all 
strearnsurfaces if the total energy i s uniform and the meridional com 
ponent of vorticity is zero upstream of the blades. Under these con- 
ditions Equation (30) becomes 



2 V 2 +f¥-oOrv - constant ( 31 ) 

In a more general case the flow may be rotational and may 
have a non-uniform energy distribution at the upstream boundary, 
the inlet. Using the subscript ( )j to denote these inlet conditions 
the ’ constants" in Equation (30) can be evaluated. 

fyf 3 - u) rv = £ , - u)(rv), (32) 

ihe functions Sj and d(rv)j are constant on each strearnsurface and 
generally each has a different value on each strearnsurface. They 
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cannot be evaluated as functions of position unless the streamsurfaecs 
are known. 

The gradient oi equation (32) is taken in order to eliminate 
the total energy gradient v £, from the equation of motion: 

7 E - 7 oOrw + 7 E, - v cd(ri/), ( 33 ) 

Since the functions and (rv)^ are constants on each streamsurface 
their gradients are normal to the streamsurface and equation ( 33 ) 
can be written as 

’£ = r U>rv + T„ -|f' - T„ (34) 

where n is distance normal to the streamsurface. 

C) Two Identities 

Two identities which can be proved for axially symmetrical 
flow using only the definition of vorticity are: 

fl • 7 tOrv = O ( 33 ) 

7 c Orv = 7 0 u)r*fl (36) 

The vector J odr has the dimension of velocity and acts 
tangentially. 

D) The Concept of Free and Bound V orticity 

The concept of free and bound vorticity provides a useful 
means of discussing the properties of rotational f iow . Although ro- 
tational flow problems are generally non-linear, they may be linear 
if the vorticity is bound. The non-linearity arises physically from 
the dependence of the solution upon the distribution of vorticity 
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which is determined, in turn, by the manner in which the vorticity 
is transported by the velocities of the solution. Some knowledge of 
the vorticity field can be obtained directly from the equations of 
motion. 

Consider first rotational flow in a region where no forces 
act on the fluid. The equation of motion is simply (from Equation 
(17b) ) : 

V*n = v-jr l/ 2 *P= vE ( 37 ) 

If the total energy is uniform on the upstream boundary it is uniform 
everywhere in the field and vZ 5 0. This is the simplest example 
of free vorticity, one in which the vorticity vector and the velocity 
vectors are parallel* . If the total energy of the fluid is not uniform, 
the velocity and vorticity vectors are not parallel. They are, how- 
ever, perpendicular to the gradient of the total energy. 

Consider next rotational flow in a region where a pseudo- 
conservative force field acts on the fluid . The equation of motion 



is (from Equation (17b) ) ■ 




V * £1 - v E- X 1 7/3 


(38) 


\V ith Equation (33) this becomes 




V */2 = v u)rw - v u>(rv), ■+ 7 E,~ Xv/3 


(39) 



♦Two parallel vectors do not necessarily have the same sense. 
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wherc 2 ^ and (rv)j describe conditions prescribed at the inlet and 
are constant on each streamsurface . Using Equation (36) this be- 
comes at once: 

( V ~ cJr) x fl = W * fl - v E, - v UHrv), - A r/3 ( 40 ) 

If inlet conditions are uniform v 2^ and v (rv) ( are zero and all of 
the vorticity is generated by the blades . In this case Equation (40) 
states that the vorticity vector, as well as the relative velocity vec- 
tor, lies in the blade surface /3 = constant; hence the term bound 
vorticity. This is the simplest example of bound vorticity. 

Forming the scalar product of the vorticity/! and Equation 
(40), the left side is zero: 

fi- r/3 = fl- j- v i E,~ u)(rv)J (4i) 

Again this shows that the vorticity vector lies in the blade surface 
if the upstream conditions are uniform. 

Consider last the flow downstream ox a system of blades. 

The equation of motion is simply Equation (40) with A = 0 : 

W * fl - p F, - v U>(rv), (42) 

If the conditions far upstream are uniform, v and 7 (rv)j are 
zero, and the vorticity vector is parallel to the relative velocity 
vector. Thus vorticity generated by a moving system of blades will 
be parallel to the Velocity measured relative to the moving system, 
regardless of the shape of the blades, and regardless of the change 
in energy effected by the blades. This is merely the "upstream" case 



-23- 



again if relative coordinates are used. 

The vorticity generated by a system of blades has certain 
geometrical properties. In the region of the blades in which the 
vorticity is generated the vorticity is 'bound' and lies in the blade 
surfaces. Downstream of the blades the vorticity generated by the 
blades is 'free' and is parallel to the relative velocity. 

E) The Transport of Vorticity 

If a fluid particle has once acquired a rotation it tends to 
maintain this rotation as it moves along the streamline. In other 
words, the vorticity is transported along the streamline with the 
velocity of the fluid. The magnitude and direction of the vcrticity 
will vary as the fluid expands or contracts, as the ve.ocity vector 

s. 

changes direction, and of course under the action of non-conserva- 
tive forces . V< ith uniform inlet conditions the equation of motion 
is given by Equation (29) with v and r(rv)j = 0. Taking the curl 

of both sides , 

PP 

~~ - V ■ r n = n • V rA * v/3 - H r- [/ _ (43) 

This shows how the time derivative of the vorticity vector depends 
on the vorticity itself as well as on the velocity and the force field . 
If the fluid is incompressible the last term of Equation (43) is zerc . 

Because of axial symmetry, only the tangential vorticity is 
associated with the radial and axial velocities, while the radial and 
axial vorticity components are associated with only the tangential 
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velocity. The tangential vorticity may be regarded as an infinite 
number of vortex rings centered on the axis. The circulation f 7 
around each ring is given by /p multiplied by the cross-sectional 
area of the ring. As the radius of the ring increases the cross- 
sectional area decreases in such a way that the mass of the fluid 
in the ring remains constant. It follows that for a given ring, f 7 
is proportional to as the ring deforms . The tangential com- 

ponent of Equation (43) expresses the law governing the time rate 



change of circulation connected with a certain mass of fluid en- 
closed in an annular vortex ring: 

< 44 > 

From the continuity equation, Equation (21): 



-L JL ru + -£L w — - PSL — 

r 3r Zz P d+ r d+ p 

With this, Equation (44) can be written: 

r p c/ 1 v_ 2 D \ _ D i p 

r P d+ 7p ~ “ * TF At * 2T At r 



/ dP 



d I 



(45) 



(46) 



The first term on the right is the rate of change in the axial direc- 
tion of the centrifugal force, which acts radially; the second term 
is the radial rate of change of the axial blade force; and the third 
term is the rate of change of the radial blade force in the axial di- 
rection. A little thought will show that these terms, in each case, 
represent moments tending to cause rotation of a fluid particle 



about a tangential axis . It is in this way that non-conservative forces 
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tend to effect a change in the circulation around a vortex ring. If 
the centrifugal force is conservative the first term is zero, and if 
the blade force is conservative A is constant, and the last two 
terms cancel. Then the circulation about the vortex ring is con- 
stant . 

Vorticity may therefore be considered as a property of the 
fluid which continually changes character as the fluid moves along 
the streamsurfaces . This concept of the transport of vorticity by 
the fluid particles gives physical significance to the non-linear 
interaction between the vorticity and the velocity and will be useful 
in understanding the various steps of the iteration process to be 
developed later. 
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V. r Ii£ V. A T HE XfA r IC A I FORMULATION OF 
THE IDEALIZED PROBLEM 

A) The Equations, the Variables, and the Boundary Con dition 

Available for the idealised problem are six equations: the 
three equations of motion. Equations (20abc); the continuity equation, 
Equation (21); the isentropy condition, Equation (18): and the geo- 
metrical relation that the relative velocity is tangent to the blade 
surface, Equation (25). These six equations relate seven depen- 
dent variables, u, v, w, p, p, A , f,--each of which is a function 
of the independent variables r and z. Clearly, since there are only 
six equations, one of these seven must be prescribed. The angular 
velocity of the rotor, cJ , is a constant parameter. 

In the direct problem the blade surface function f(r,z) is 
prescribed and there are six equations for the remaining six de- 
pendent variables . The necessary and sufficient boundary conditions 
are the same as for flow through the same channel with no blades act- 
ing except that a sort of Kutta condition must be applied at the trail- 
ing surface of the blade region. The boundary conditions will be 
discussed in more detail after the number of dependent variables has 
been reduced . 

In an indirect problem where the local rate of energy input, 

A (r,z), is prescribed there are again six equations and six depen- 
dent variables. The question of blade existence, so important in the 
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indirect problem with force components prescribed, does not arise 
here, for with the force field expressed in terms of the two functions, 
A and f, blade existence is assured. The boundary conditions, how- 
ever, are less understood than in the direct problem, particularly 
the conditions for f(r,z). It is to be noted also that while the direct 
problem, with f prescribed, can be formulated as a problem in the 
meridional plane independent of A , the indirect problem, with A 
prescribed, is not independent of f . The two problems are therefore 
quite distinct. 

The direct problem is formulated below in terms of a stream- 
function for the meridional velocity. The discussion of characteris- 
tics and general properties in Part VI also refers to the direct prob- 
lem . 

B) The Differential Fquations for th e Streamfunction 

The continuity equation is identically satisfied by the stream- 



function ip , defined by: 



u _ _£~ j_ 

a* p r oz. 

IV P* J_ 0 'P 

a*. ~ ~ r r Or 



( 47 ) 

( 48 ) 



For cases in which the fluid is incompressible take 1 ] —£•* = | . 

The equation of motion is 

V * fl = 7 E -A r/3 ( 38 ) 



=# In this definition a* and p* are considered as constants which de- 
scribe the initial total energy. If the initial total energy distribution 

is non-uniform, a* and p* are constants which describe the total 
energy at some representative point cn the upstream boundary. 
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Up3trearn or downstream of the blades the force magnitude A is 
aero and since the total energy remains constant on each stream- 
surface 

17 £ = ^ M ( 4(? ) 

Written in intrinsic coordinates Equation (38) becomes 



_ T OE 

~ ,f7 On 



O V q s 

-fi-n n. ^-5 

Using the definition of <Tl n and fl s , Equation (3), the normal 
component is 

na - -id •LS y * E 
'l r dn ~5n 

The other components show that is zero . 

o S 



(30) 



(31) 



.'ince the total energy is constant on each streamsurface , 
Equation (34) can be replaced by 

7£ ~ 7n 4n “ ) [( r ' / k ~(rv),] + T n ~ ( 5 ?.) 

and Equation (31) can be written for the upstream and downstream 
regions respectively as: 

n q s .-L/rv),^h.-.^l (53) 



q<=7s= T* [(r\z) t -^r 2 ] 



0 ( rv)t 



+ dE, 



(34) 



On dn On 

In the region of the blades \ is not zero and the total energy is 

not constant on each streamsurface. *n expression for the tangential 

vorticity can be obtained from Equation (41): 

n • ^ = n • f T n ±.[e, - unrv),j (41) 



Expressing the left side in cylindrical coordinates and the right side 
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in intrinsic coordinates there results : 

-Sfr &[ (55) 

But fl n = - ~pr jlf' , Equation (3); and X = q s , Sq. (22b) 
so that 



n ' 0 CX 

*\n - r ZEU — 



X 



<=fs 



a ry 

os 



rq s 



(56) 



and Equation (55) becomes 

7 = K & - /> &) - 4 & [e, - uxrvl] (57) 

Equation (5 7) expresses the tangential vorticity Q in the region of 

the blades. It is significant that this relation is independent of the 
blade force magnitude A . 

All of the three basic equations, Equations (53), (54), and (57), 
are non-linear if the upstream conditions, E^ and u)(rv) , are non- 
uniform, or if the fluid is compressible. If the fluid is incompres- 
sible and and (rv)j are constant, the equation upstream is linear. 
It is simply i\= 0, but the downstream equation, Equation (54), is 
still non-linear. However, the equation in the region of the blade. 
Equation (57), actually is linear if the fluid is incompressible. For, 
from Equation (25) 

rv = cOr z v- r 2 (-f r u (58) 

so that in the region of the blades 

7 =/£ fr - f r§z){ + r z {f r u * £ 

<7s u J 



( 59 ) 
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Ncte that the operator (^zj- r ~ represents differentiation 

along the trace of a blade in the meridional plane. Equation (59) 
is clearly a linear equation for u and w when the last term is Eero 
and f(r,z) is prescribed. 

If q is written in terms of the 3treamfunction in Equations 
(5 3), (54), and (59), the three differential equations for the stream- 
function are obtained. By direct substitution, using Equations (47) 



and (43) , 

a u o w 
'l ~ 02. ~ or 

_ fi. 

P 



. a r 

f / O'P ; 0* 'P _ d±J_OP d_ f_J_ bp 7 

> r [ <Dr z ~ r Or + dz z or p or 7 oz P Oz J 



The last two terms can be written more simply after a change of 



independent variable as 

0±_J_Of_ 0^ ! BP _ 

Or P Or Dz P Oz. ~ on P On 

The derivatives taken normal to the streamsurfaces of 



(61) 



functions which are constant on the streamsurfaces can be written 
in a very useful form by considering the streamfunction as an 
independent variable. For instance, 3ince rv is constant on each 
streamsur face downstream of the blades , (rv)^ depends only on the 

streamfunction tp . Therefore 

i>(rv) t _ djrvk O'P _ c/(rvk / Qs P„) (£> 2 ) 

Dn d<p On d l a, fi* J 

and similarly for any function which is constant on the streamsurfaces 
Using the relations developed in Equations (6l) and (62) and 
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the definition of <p the differential equations for the streamfunction 
can now be written for the three regions of the flow, equations (53), 
(54), and (59) then become: 

Upstream of the blades: 



+ = { Pn( lnF )n 




2 cJ £, 

dp a* 




( 63 ) 



In the region of the blades: 

( t+Z 2 ) 'f’rr ~2RZ ip ri + (l+R 2 )ifi zz + (ZZ r -f?Z z -pJ b '(ZR r -RR 2 )^J z 
= RZ4*z^(inf) r +[(l4R*)') l z.-RZ i> r ] (lnf) z 

+ /jLf r z°L 7 ( 64 ) 

( rj L a? a m * J 

where R = rf^ , Z - rf ^ , and the subscripts r and s denote partial 
differentiation . 



Downstream of the blades: 



^ rr ~ -p <P r *■ 4>z. z - 




(65) 



The equations are rewritten below for the flow of an incompressible 
fluid with uniform inlet conditions in order to discuss the rotation- 
ality effects separately from the compressibility effects. For the 

Z? 

incompressible fluid, a* and are taken as unity. 

Upstream of the blades: 



'Prr - 7 * ° 



(66) 



In the region of the blades: 

(i+Z*) *rr~ 2 RZ ifj rz +(/ + R 2 ) <P ZZ 



+ (ZZ r -RZ z - -f)^ -(ZR r-RRz)^ * 2 <ArZ 



(67) 
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Downstre am of the blades: 

trr -y'Pr +4>zz= - [l^h ~ o)r2 J^ (r '' > * ( 68 ) 

With uniform inlet conditions the flow upstream of the blades 
is irrotational . equation (66), which is clearly a linear partial dif- 
ferential equation, states simply that the tangential vorticity n is 
zero. In the region of the blades and downstream of the blades the 
flow is rotational. However, in the region of the blades the fluid 
is constrained to movement on the blade surfaces , and the vorticity 
is 'bound' and lies in the blade surfaces. Because of this constraint 
the equation for the streamfunction in the region of the blades, Equa- 
tion (67), is linear when the blade shape functions, R and Z, are pre- 
scribed. The equation for the downstream region, Equation (68), is 
a non-linear equation for ^ because ^ appears as an independent 
variable on the right side. However if the right side is a known 
function of position the equation is linear , although non-homogeneous , 
and can be handled in the same way as Equations (66) and (67). This 
suggests an iteration process in which the right side is estimated and 
the equations are solved. A new right side is calculated from this 
solution and the process repeated until the desired accuracy is ob- 
tained. Actually (rv)j is a function of only and the non-homogeneous 
part, the right side of Equation (68), can be evaluated as a function of 
ij' and r for all the downstream region once the flow conditions at the 
blade trailing edge are determined. There are several reasons why 
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this iteration process is a practical and relatively simple method 
of solving this problem. I'he nen-homogeneous part of Equation 
(68) has a simple physical meaning. The right side is simply rq , 
the moment of the tangential vorticity about the axis , and is pro- 
portional to r p . It was shown that under the action of a conserva- 
tive force field the circulation p remains constant along the stream- 
surfaces (Equation (46)). This suggests that a gooJ first approxi- 
mation might be to assume that the circulation is constant on the 
streamsurfacea of an irrotational flow through the same channel, 
thus neglecting the interaction of q and the other vorticity compon- 
ents . actually a much better approximation can be made by evalu- 
ating the function (rv) as a function of ij; a * the blade trailing edge 
using the blade geometry and the meridional velocity of an irrota- 
tional flow through the same channel. Once (rv)^ is a known func- 
tion of the right side is a known function of ip and r. Here the 
interaction of the vorticity is accounted for, but the blade trailing 
edge conditions, which define rv(»jj ), are of course approximate. 

In the numerical method of finite differences used in the following 
examples it is only necessary to 3olve completely the final step of 
the iteration process. Each step only need be carried far en„_gh 
to assure an improvement of the estimate of the right side . The 
application is relatively easy when the relaxation technique is used 
since the exact status of the solution is clear at all times. 
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Ahen the fluid is compressible all of the equations are 
non-linear. The term on the right side, (in > represents 
the entire effect of compressibility if the flow is irrotational . For 
rotational flow another effect of compressibility is to modify the 
effects of the rotation. Thus the density ratio ) is a multiplier 
of the rotationality terms in liquations (63), (64), and (65). Again 
an iteration process is used, in which the term * s est i- 

mated and a solution is obtained. This solution is then the basis 
for a new estimate of the density term, and the process is repeated 
until the desired accuracy is obtained. * 

The conditions under wuich the above iteration processes 
are convergent have not been rigorously established. For flow in 
which a certain governing velocity is supersonic convergence is 
questionable. If stagnation points occur in the fluid the question of 
convergence is connected with the question of proper boundary con- 
ditions and the fact that the streamsurfaces are characteristic sur- 
faces of the flow. Convergence will be discussed separately for 
each example solution. 

The equation for the flow in the region of the blades is 
somewhat simpler than shown in Equation (64) if radial blades are 
prescribed. Radial blades, blades which are generated by radial 
lines passing through the axis of rotation, are necessary in very 
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high speed machines because of the high centrifugal forces. The 
trace in the meridional plane of a radial blade is of course a radial 
line, hence f(r,z) is actually f(z) only. Then for radial blades: 

R = O, Z ~ r (z) , Z r = f z (z)j e-Z-c.j (69) 

and Equation (64) becomes: 

(l + Z’j+rr-ll-Z^-k'l'r 

+(l*Z’l'l>,(lndr + 'h(l><e)z +/£)’ r *a?[4fr- !i 5$ < 70 > 

The significance of the simplification is that the cross derivative, 
if) > does not appear when the blades are radial. 

The most general form of the differential equations for the 
streamfunction is that shown in Equations (63), (64), and (65). These 
are the equations for the isentropic, rotational flow of a compres- 
sible fluid acted upon by an infinite number of arbitrary blades. 

The flow may be rotational and the energy distribution may be non- 
uniform at the upstream boundary. An energy relation which ex- 
presses density in terms of the gradient of the streamfunction and 
the total energy of the fluid must be considered with these if the fluid 
is compressible. 

C) The Isentropic Energy Equation--Density as a Function of £/ass Flow 
The energy equation, Equation (32), can be written as 

V 2 -+- -j 2 j- = cOrv ■*■£,- od(rv), (71) 

Let a Q be the velocity of sound when the velocity V is zero. The 
constant a Q will have a different value on each streamsurface if the 



energy on the upstream boundary is non-uniform. 



.. non-uniform 



v.ill be denoted by a prime, e.g. a ' . Then equation (71) becomes 



a z / *-' ( V - (rv- rv,) ) 

a' 0 * ~ ' 2 l ai* J 



(7s) 



or in terms of density 

P _ ft y-/ / u z . 1 r-/ ( vj_ tArv - <£rv, \ 1 Y-> 

Po ~ j 2 '&<!* G-'S) 2 l oq 2 *• «o l )J ' 

Introducing the streamfunction defined by Equations (47) and (4S) , 



^ _ 



P' 

f o 
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his equation gives the density in terms of the local mass flow 

fa . / P* 4 >f+ Pi 



2 2 

Pqs 

P' 2 a ‘ 2 



ct' z Pc'* 



and the effective tangential velocity 
Veff. _ v 2 - ZiA)lr/-r\f,) 



(75) 



„ , (7b) 

I a l 

This last term includes the effect of changes of total energy ac a re- 
sult of blade action as well as the actual effect of the tangential ve- 
locity . 

The solution of Equation (74) is presented graphically as com- 
putational curves, Figure 42. 

!’wo other useful forms of the energy equation can be derived 
from Equation (72): 

* ■ fftfJll ‘ 2 

P ! P:J r 



Pc 2 



df_ _ 
P 



2 U)(rv- rv,)- v 2 



-*-f&! .£l 

y-/ f a* a.*J r 

' ll'Pr+'Pz*'\ ' ,/ 2. a 2 2 o)(rr- rr,)-v* j 

2 d l T 2 1 P&IY-/ aS ct * J 






2 f Op Off) 2^(rv-r 

r-i I aS ~ a '*! 



'-rv f )-v 2 



(77) 



- (70) 
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D) The > etching of olutions in Adjoining Regions 

There are three distinct regions of the flow corresponding 
to the three differential equations, Equations (63), (64), and (65), 
and solutions in adjoining regions must be matched on the connect- 
ing bounuaries. The streamfunction itself is always continuous 
across both the leading edge and the trailing eage boundaries. The 
I'utta condition imposed at the blade trailing edges to make the so- 
lution unique requires that the velocity and pressure be continuous 
there. Consequently, the matching conditions at the trailing edge^ 
are that the streamfunction and its first derivatives, as well as the 
tangential velocity, are continuous across the trailing edges. This 
provides the means of evaluating the function (rv) t for the down- 
stream region. Liihce rv is continuous it may be evaluated just up- 
stream of the trailing edge by Equation (25), which can be written as 
rv = u) r 2 + r 2 {f,.u + (79) 

The matching conditions at the leading edges are generally 
much more complicated. In the real problem with a finite number 
of blades the fluid flows smoothly through the blade system, .i 
"stagnation point" occurs at some point near the leading edge and 
the fluid flows smoothly off the trailing edge. The flow is similar 
to the flow about an airfoil. However, in the idealized problem the 
number of blades is infinite, and the velocity at the leading edge 
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will generally be discontinuous. The flow is vary similar to the flow 
through an infinitely closely spaced lattice of airfoils, and the lead- 
ing edge discontinuity is the same as that occurring at the leading 
edge of the lattice. In problems of practical interest the pressure 
and tangential velocity will generally be discontinuous at the leading- 
edge . Continuity of the first and higher derivatives of the stream- 
function depends on the shape of the blade near the leading edge. It. 
can be shown that the first derivativesof the streamfunction arc edn- 
tinucus across the leading edges provided only that the leading edges 
lie in meridional planes (Figure 4). 




It i3 supposed that the tangential velocity v jumps diocontinuously 
across the leading edge. This implies that a discontinuity in density, 
pressure, and velocity may also occur there. The first derivatives 
of the streamfunction are proportional to pq $ , the local mass- velocity 
in the directions of the derivatives . pq^is of course continuous on 
both sides of the leading edge, and continuity of mass require that 



-44- 



(P,)(ch), =(P*)M* in Figure 4. Therefore pq^is continuous every- 
where, and the first derivatives of the streamfunction are continuous 
everywhere. Tor blade surfaces which satisfy this condition, the 
matching conditions at the leading edge are that the streamfunction 
and its first derivatives are continuous there. 'The tangential, ra- 
dial, and axial velocity components and the pressure and density 
may be discontinuous. If the blades do not satisfy this condition 
the streamfunction derivatives will generally be discontinuous. 

3) The Boundary Conditions 

-» typical region of flow is shown in Figure 5 . The axially 



Do iv nsf r* a m 




Figure 5 

A Typical Region of Flow 

symmetrical surfaces AB and CD are boundaries of the channel and 

» 

are streamlines of the flow. The upstream and downstream boun- 
daries are the surfaces i:C and BD respectively. The region of the 
blades or the body force field, over which the blade surface function 
f(r,z) is prescribed, is shaded. 
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I? or each atop of the iteration the equations for the stream- 
function, Equations (63), (64), (63), are non- homogeneous linear 
second order partial differential equations. The necessary a:iu suf- 
ficient boundary conditions for this problem are well known. I he 
streamfunction </j , or a linear combination of its derivatives , 
and p^ , must be prescribed at every point on the boundary. The 
boundary conditions on the streamfunction for a typical region of 
flow (Tigure 5) are: 

p = constant on surface AB 
p - constant on surface CD 
P = p (r) on surface 
= 0 on surface BD 

The functions v, p, and p, which can be prescribed ^nly once on each 
strearnsurface , are prescribed at the inlet, on surface - . 

it the leading edge the matching conditions arc that the 
streamfunction and its first derivatives are continuous. the 

Cr 

trailing edge the s t re amf unction and its first derivatives, as well 
as the pressure, density and tangential velocity, are continuous . 

This is the ivutta condition. 

The iteration process is based on an elliptic partial differ- 
ential equation. In ~^art VI it is also shown that the complete non- 
linear differential equation, which includes the corr.pressioility 
term, becomes 'hyperbolic” when tne meridional velocity or, if 



blades are present, the relative velocity, becomes supersonic. 

The boundary values for this case are quite complicated and solu- 
tions may not be unique. In all the examples solved here the gov- 
erning velocity is subsonic, although the total velocity may be 
transonic, as in Part X. 

F) The Limiting Flow F ar Downstream of the Blades 

If at a finite distance downstream of the blades the channel 
boundaries become concentric cylindrical surfaces and if these 
surfaces extend unchanged downstream to infinity, general quali- 
tative statements can be made concerning the limiting flow at in- 
finity. It is supposed that the total energy of the fluid is uniform 
and that the vorticity and tangential velocity are zero on the up- 
stream (inlet) boundary. The energy of the fluid is changed and 
vorticity is created by the action of an infinite system of blades 
rotating at a constant angular velocity , Nothing is said about 
the channel boundaries upstream of the cylindrical region or about 
the shape of the blades. 

It is clear that in such a region the flow is independent of 
the axial coordinate and the radial velocity is zero. Then radial 
equilibrium of the centrifugal force and the pressure gradient 
force requires that (from equation (20a)) : 
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The vorticity, being generated by rotating blades, is parallel to 
the velocity relative to the blades. Therefore from Equation (42): 



d I'V z _ o u) r - V dry 
dr r dr 



(81) 



From Equations (80) and (81) and the isentropy condition, Equation 
(19), the following relations are derived: 



/ c/P J_ V z 

f dr p d z 



(82a) 



<d_a_ s 

dr 




(82b) 



2 

d M* _ £ _v_ 2 / , cor d ry \ 

dr r a* l dr j 



(82c) 



d M 2 _ <£? V z / ! t-l m z dfv\ 

~dr ~ r a *{ ‘ ^ 2 ' ' ~ dr J 



The energy equation, 

/ + - 



Equation (72), can be written as 
2 

&o , y-/ t*>rr 

a 2 - 2 a z 



Using this , 



d M 2 _ 2 V 2 *"'/ / JW 2j£r_s_j 

dr ~ r a 3 - a 2 - i 2. af ) 



(82d) 



(82e) 



If rv is a known function of r or of in this region the flow 
can be completely determined by the above equations. It can be 
seen directly from aquations (80), (82a), and (82b) that the pressure, 
the density, and the local speed of sound are minimum on the inner 
boundary and increase monotonically to a maximum on the outer 
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boundary. This general result is completely independent of the blade 
shape -nd the channel shape, it the blades are stationary it follows 
from equations (82c) and (32c) that the total velocity and the corres- 
ponding ucu number are muxinnum on the inner boundary and de- 



velocity may have a maximum or minimum anywhere, depending on 
rv, but if rv is monotonic then the axial velocity is also monotonic. 

In order to examine the total mas 3 flow the channel is given 
the same cross-sectional area far upstream and far downstream. 
The total energy of the fluid is constant. If the tangential velocity 
is zero upstream and constant do w nstream the total mass flow equa- 
tion is 



where ( ) and ( ). denote values far upstream and far downstream 

a L 

respectively. 

Using the isentropic relation. 



crease monotonically to a minimum on the cuter boundary. The axial 
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and integrals of equation (82), Equation (83) becomes 




2 



rd r (84) 



where 




#7* is the maximum which occurs at the inner 



radius , taken here as unity. 



.he solution of 'equation (04) is shown in figure 43 where 
















I 
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2 2 V ^ 

/ J is plotted as a function of Mf for several values of r . 

1 < 377 *; 

These results are expressed in terms of local Mach numbers in f igure 
44. where the axial Mach number /-£-) 2 , is plotted as a function of 

V & ' Zf 



1/ Z 

the inlet Mach number M, for several values of 



. It is seen 



that for a given inlet Mach number, M, 2 » and a given deflection, 

there may be two solutions for the maximum axial Mach number 

downstream, but if P~) is above a certain critical value there will 

( a *- 'z, 

be no solution. Vhen two solutions exist one value of the axial ve- 
locity may be subsonic and the other supersonic, or both values 
may be supersonic. It must be understood that ('a~) x/ - an< * (£1 
are maximum values on the downstream boundary and that when 
these are supersonic the minimum velocities, and, even the mean 



velocities, may still be subsonic. 

These results depend only on the limiting flows far upstream 

$ 

and far downstream and are independent of the intermediate flow. 

It is assumed, however, that the intermediate channel and the blades 
are such that an lsentropic transition from the upstream flov.’ to the 
downstream flow is possible. Nevertheless, this analysis shows 
that for a given inlet Mach number there is a maximum tangential 
velocity that may be generated by stationary blades. No isentropic 
solution exists for tangential velocities exceeding this maximum.. 
Similar results would be obtained for distributions of tangential 



velocity other than constant, and for rotating blades. 






I 



-50- 



VI. PROPER FIEd OF TH C TRANSONIC FLO / 
OF A PERFECT FLUID UNDER THE ACTION 
OF PSEUDO-CONSERVATIVE FORCE FIELD 



a ) The Characteristics of the Problem-- The Governi ng Velocity 
a he applicable equations upstream and downstream of the 
force field are shown in Equations ( 19 ), (20), and (21) with A = 0. 
Adding to these the four differentials for u, rv, w , In p , and using 
Equation (19) to eliminate the pressure terms, eight algebraic equa- 
tions for the eight partial derivatives are obtained: 




r 




O 




c/ 

r 



(S3) 



Thj diiferential equations for the characteristics are obtained when 



the determinant of the coefficients is 3 et equal to zero: 






( 
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u 

0 

0 

1 

dr 

0 

0 

0 



v; 

0 

0 

0 

dz 

0 

0 

0 



0 

a 

0 

0 

0 

dr 

0 

0 



0 

w 

0 

0 

0 

dz 

0 

0 



0 

0 

u 

0 

0 

0 

dr 

0 



0 

0 

w 

1 

0 

0 

dz 

6 



u 

cl 

0 

0 

u 

0 

0 

0 

dr 



0 

0 

a 2 

w 

0 

0 

0 

dz 



= 0 (86) 



Expanded, this equation is: 

w z {w 2 a 2 )dr*+ 2 uw (2 a £ )o/r 3 dz. [6 u 2 w z -(u 2 + w 2 ) <x s ] d r* dz* 

+ 2 u.w [2 u z - a*) drdz 3 + u z ( u L - cf) d r 4 = O 

which on factoring yields: 

< 87 > 

dr uw ± ajct*r w* - a"*- 

~dl = (OS) 

-Equations (87) and (68) are the differential equations for the physical 
characteristics of the problem. These equations represent, respec- 
tively, the streamsurfaces and the 'meridional Iv.ach surfaces". 

Thus the character of the flow upstream and dosvnstream of the blades 
depends on the meridional velocity, not the total velocity. That is, 
the problem is "elliptic' with respect to the variables u and w if the 
meridional velocity is subsonic, and is "hyperbolic" if it is supersonic. 
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The velocity which determines the character of the problem is 
called the "governing velocity", thus the governing velocity up- 
stream and downstream of the blades is the meridional velocity . 

In the region of the blades the applicable equations are 
those given above but with A f 0 . In addition Equation (25), which 
states that the relative velocity is tangent to the blade surfaces, 
must be included. Equation (19) can again be used to eliminate 
the pressure and Equation (20b) can be used to eliminate A . 
Because of the constraint imposed by the blade surfaces (Equation 
(25)) it is also possible to eliminate the angular momentum rv. 
The blade surface function, f(r,z) is prescribed. Then, including 
the differentials for u, w, and In p, there are six linear algebraic 
equations for the six partial derivatives: 





U_ 

r 




du 



(39) 



d w 




where 



G = £ UH-u + u Z £- rR 



4- U\/V 





rc tic 
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.ctting the aeter r.-iino nl oi the coeflieientsequ&l to z. : 
following equation is obltineU: 

wf(/ + r*+z*)w 2 - 0 + R v a z j dr 3 - {3(hR*+2 z )uw* [(UR 2 ) a -2 RZ^Ja z } dr *dz 
+ [3(Uf?*-*Z 2 ) -[(/ +2 7 )w -2 RZ uj a z }dr dz z 

- uf(/+R 1 +Z 3 )u 2 -f/+z z )a z }dz 3 -0 (90) 



On factoring this yields: 

W d r - u dz = O , and 

f/J^ Z Z w a a *)(,+/?*) dr 2 - UW + ) JhR ~ 3 Juz* drdz 

I !+«*■ I WTFJITz* Jur* JiZz 3 ) 

+ (i±zh£ U g - a*j(i+Z z ) dz z - O 

The differential equations t>r the physical ch-v: «*cfcori^tics are 

therefore: 

u 



c/r 

dz 



w 



( 91 ) 



dr _ fUR^*Z^jjddZ_ +&Za*± gJfUR^z^ia 3 ** 3 * (RutZw) 2 - a*J c ? 



dz 



(/+ R*+ z*J iv *- (/ + R S ) a 2 



( 92 ) 



where Ru+Zw is oh a relative tangential velocity, (v~u>r). 

The characterisrics oi the flow in the region of the blades 
are the streamsurfaceo ( .quation 91) ana tne surfaces given by loca- 
tion (9-2). /heth^r or not the char acteriatics given by Equation (9d) 
are real depends on the total rel^-ive velocity, not on the meridional 
component of the velocity cts in the upstream and downstream regions, 
’he governing velocity in the region of the blades is therefore the 
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VGlocity relative t.., the blades. 

The form of ^auction (9*2) suggests the following transform- 



ation: 



U 



J77z* 




J7^r^z^_ 

JITr* 



uv 



dr - JT+r* dr ; dz = /TTz 2 - dz 



(93) 



After this transformation Equations (91) and (92) become: 



dr _ J± 

dz w 



dr_ 

dz 



UW + 



2#Za^_ 

JTtH 2 - JI+z z 



±0/5 



2 _ 2 
tw + 



(,+pvo+z 2 ) u 



- a 



a 



(9i) 



This is no great simplification except when R or Z is zero ever/where. 
For radial blade R is sere, the transformation is: 

Z7 - u ; w - / / TiE* w 



dr - dr ; 



dl = 



(95) 



and the equations for th>- caaracteristics , in terms of the transformed 



variables are: 



dr _ tu 
dz “ "S' 

dr _ uW±ctJ CL** ifr 2 - Cf 2 
dz. W 2 - ot 2 - 



(96) 



Thus the characteristics for the transformed problem are the same 



as for two-dimensional flow with velocities TI and w in the r, z space. 



h>) The Complete Kon- linear Forms of the Equation 



With uniform inlet conditions the equation for the gtrean.- 
function for compressible flow downstream of a system of blades 



I 



I 



I 



relating at angular velocity tO i» (from liquation (65)): 




( 97 ) 



The density derivatives can be written using Equation (78) as 





Then, eliminating the density derivatives from Equation (97) and 
introducing the velocity components, the complete non-linear form 
of the equation for the downstream region is obtained. 



The noil-homogeneous part, "the right side* , represents the effect 



connected w.th the fact that part of the total energy occurs as kinetic 
energy due tc the tangential velocity. It is clear that the effect of 
rotationality becomes aero if the meridional velocity is sonic, and 
that the rotationality effect on one side of the sonic line is exactly 
opposite the effect on the other side. 

The equation for the streamfunction for compressible flow 
in a region of radial blades is (from Equation (70)): 




w 2 

of rotationality. The terms containing the derivatives of — are 
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The density derivatives, from Equation (78), can be written as: 



/ DP 



co 
a ; 






/ ~ 



u^+O+Z*) yy z 



J_ £f _ jJzgjifat' 

p 3Z 



<1^1) 



/_ (H-Z x )w* 

• 9 



The complete non-linear equation for the region of the blades is ob- 



tained by eliminating the density derivatives from Equation (100) and 
introducing the velocity components in the coefficients: 

(/-■^)(/*Z , )4’rr + ^£f'+Z I )<Prz + [l-ll + Z 2 ^] <Pzz 

-{(l-Z*)+ IUZ 2 )~r^[2-^^ U+Z')„']rZ&} ± 



+{-£z§r) = 2--^ Z[l- ] ( 102 ) 

C) The Cushioning “ ction between Subsonic and Supersonic Regions 
The complete non-linear equations (Equations (99) and (102)) 
are obviously quite complicated. It is believed, however, that the 
transformation suggested (Equation (95) ) would lead to some simpli- 
fication of Equation (102). Unfortunately, time has not permitted 
a more thorough investigation of the possibilities of this transform- 
ation . 

Examination of Equations (99) and (102) and the characteris- 
tic equations, Equations (88) and (92), reveals that the non-homo- 
geneous part, the right side of the equation, changes sign whenever 
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the character of the flow changes from "elliptic" to "hyperbolic" . 
Furthermore, for flows of practical interest the vorticity term 



(f_P \ Z [ £1 _ — 

I* [ a+ &+ J c/tp 



£? L , or Z — 4- Z 

a* t J a* p* 



will have the same sign throughout the flow region. The right side 
of the equation can be regarded as a forcing function acting to dis- 
place the streamlines. In regions where the flow is entirely sub- 
sonic or entirely supersonic the forcing function will have one sign 
throughout the region. But if the flow is transonic, i.e. part sub- 
sonic and part supersonic, the forcing function has one sign in the 
subsonic part and the opposite sign in the supersonic part. There- 
fore the deflection of the streamlines, brought about by vorticity , 
will be less if the governing velocity is transonic than if it is en- 
tirely subsonic or entirely supersonic. Consequently the meridi- 
onal flow will be "smoother", for the vorticity in one region tends 
to counteract that of the other region. It i3 believed that this mu- 
tual cushioning effect is the explanation of the phenomenal effic- 
iences observed in compressors in which the relative velocities 
in a region near the tip are supersonic. 

The fact that the forcing function, the non-homogeneous 

part of the equation, contains a factor of the form (l-Iv *') , where 

0 

Ivig is the "governing" ta&ch number, leads to the following general 



conclusions: 
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1 . The deflection of the streamsurfaces induced by a given 
strength of vorticity at a certain point in the flow is aero 
when the governing velocity is sonic at this point, has one 
sense when the velocity is subsonic, and the opposite sense 
when it is supersonic. 

2 . The deflection of the streamlines brought about by vor- 
ticity in a region is less when the governing velocity is 
transonic in the region than if the velocity is entirely sub- 
sonic or entirely supersonic. 

These conclusions are confirmed in all of the examples of 
compressible flow which follow. In fact, it seems that even stronger 
statements should be made. A comparison of the deflection of the 
stream9urfaces brought about by vorticity indicates that the deflec- 
tions increase from aero to a maximum and then decrease as the 
vorticity (the strength of the actuator) is increased from zero. 

The governing velocity is subsonic in both examples . 
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vil. fhE . TION OF TH'C FINIT - 

ERdNCE PRO BIGN 

The finite difference problem Gitfers from the differential 
equation problem in two ways: in the finite difference problem, 
the boundary conditions are prescribed at a finite number of points 
on the boundary and the desired function is to be determined at a 
finite number of points within the boundary; the differential equa- 
tion is essentially replaced by a number of simultaneous approxi- 
mate difference equations which are themselves only approximately 
solved. However, the set of points at which the desired function is 
sought can be made very dense, and the system of difference equa- 
tions can be solved as accurately as desired. Therefore, barring 
unaccounted for singularities, any degree of agreement between 
the two solutions may be obtained. The two problems are neverthe- 
less quite distinct. The finite difference problem is solved by outh- 
wall's relaxation method (Reference 8), which is defined as "a sys- 
tematic sequence of localised changes of the wanted function that 
steadily brings the 'residuals' toward their desired value.” 

The fact that the physical aspects of the problem are always 
evident and that the status of the solution is always apparent is a 
great advantage of the relaxation method. This is a great aid in 
solving the non-linear turbomachine problem, for the action of the 
vorticity and the compressibility is apparent as a force" which 
causes the streamlines to deform from their incompressible , irro- 
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tational positions . 

The power o i the relaxation method i.-., indicated in the fol- 
lowing quotation from Southwell (Reference :•>, p. 3): 

' (The) use of finite differences is no. new, nor (is the) 
evaluation of the wanted function at nod<.d points of a regu- 
lar net. but in concentrating attention on the data, an- in 
recognising that these are never exact, they subordinate 
mathematical to physical aspects in a way that can alter 
drastically the course of a theoretical research. Z’is- 
carding orthodox for relaxation methods, an investigator 
finds his outlook quite transformed: full scope remains 
for ingenuity and special artifice, but any problem that 
can be formulated can be solved.’ 

It is not clear how solutions which are unstable can always be found 
by the method of finite differences. Intuitively it seems that if the 
configuration is unstable the disturbances inherent in the relaxation 
process may be amplified by the relaxation, thus making converg- 
ence questionable . 

A) The Difference equations 

The most general form of the differential equation for the 
streamfunction is that shewn by Equation (64): 

A ip rt - f 2 B i^ rz +■ C Lpzz +0 = 0 

HO 3) 
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where /• ,B,C are functions of r and z, and D is a function of 
\p r , ip z , r, and z. As written for e«ch step of the iteration v itL 
the non-nomogeneous part known, Equation (64) is an elliptic equa- 
tion, for 

AC - B 42 = l + R 2 +Z 2 >' O , 



and can therefore always be transformed into the normal form 

(ft xx + 'Pyy 4 ^ * j R j Pr j <Pz I — C> (- 104 ) 

,Vhen the cross derivative term ip r2 occurs, as in Equation (64), the 
difference equation may be written for either the original equation 
or the transformed normal equation. Because of the added difficulty 
of matching solutions when one region has been transformed, the 
difference equations are written here for the original form. 

If the screamfunction p , which is a solution to a partial 
differential equation such as Equation (64), is regular in some neigh- 
borhood of the point (r , z ) , it may be represented in this neigh- 

tj O' 

borhood by the Taylor’s expansion: 

P'(CZ) = ^P 0 + <p r (r-r 0 \+ cp z (z-z 0 ) 

4 [Prr (r-r 0 ) Z -f- 2. ip rz (r-ra)(z-z 0 J h tfi Z2 (z-z 0 \ Z ] 

^ ' (165) 



where <p 0 is the value cf p at (r Q 

evaluated at (r , z ). The re - 
v o o' 

fore if (r 0 ,z^) is one point cf 
square lattice with an in- 



i. ) and all the derivatives are 
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terval between points of S , 
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the value of p at the surrounding points can be expressed by moans 
of equations (105). If fourth and higher powers of S are neglected 
the result is a system of simultaneous linear algebraic equations 
which cun be solved for the derivatives to give: 



[p r is. 

</Va « 

4>n « 



«= ^ Rz - 

S* / Pzz 

ts-Pt+Pr- 

4 S 2 



2 S 

p + p3 -_2_h 

S z 



( 106 ) 



Cn eliminating the first order derivatives it follows that: 

5^/ * 02 * 9^.3 + ~ 4 f 'Prr + Pzz )S Z 

Ps * Pi * P7 * Ps - * Po - (Prr-i- Pzz) S * 

The difference equation is obtained when Equation (106) is 
substituted into the differential equation. Thus for Equation (64), 
which is 
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/ /+Z 2 ) <jj tt ~2f?Z p rz +( '/+R 2 )p zz + fZZ r -RZ z -j) tp r - (ZR r -RR 2 ) p z = Q ( 1 08 ) 

the difference equation is: 

[(hR^-tiRrRRji-] * [(i*z*I~(zz,-rz ! -t)t] 4>a 

+[0+R , }'(Zfi r RR l )&-] ,Jj 3 + [f/*Z 2 J- fzZr-PZi-i}^ J #4. 

- -t-RZ(p s -p^ *p 7 - i/j 3 j - 4-[ /-^ j- (R*+Z 2 )] 4j q = Qs 2 ' (109) 

where 

Q'ZZTtIz *[('+Z’l4’r-*Z'l>.](lnP) r + [(/+R’)4‘ 1 -l?Z4, r ](lnp)z 



-P ) 2 r 2 d f tdrS, I 

r dpi^~af) 
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,.hen the blade shape function f(r , 7 ) is prescribed all the coeffic- 
ients can be determined. Equation (109) i9 then the difference equa- 

2 

tion applicable in the region of the blades. The residual, <rj_ S , is 
a complicated function of position which is estimated for each step 
of the iteration by means of the solution of the previous step. 

In the same way the difference equation upstream of the 
blades is, from Equation (63): 

4>, +(l'h)'Pz + <P 3 ( U °) 

where 



* • *" '[' m £ 

The difference equation for the downstream region is the 
same as Equation (110), but the residual, from Equation (65), is 



given by 



= 




( 111 ) 



The difference equation for each of the three regions is 
shown in Equations (109). (HO), and (111). A general form might be 
written as 

DiM = S 2 <% (112) 

Here Dj(<t>) represents the left side of one of the difference equa- 
tions written for the i^ 1 point of the net and Stf. is the residual 
there. Dj( ) is the "difference operator' , and is itself a function 
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of position. That is. the coefficients of the left side of the differ- 
ence equations have different values at each net point and these 
values depend on the prescribed shape of the blades. 

If the blades are prescribed as radial, the expression for 

D- is more simple than in the general form (Equation (109) ). For 
radial blades 

D /W = +[f!+Z 2 ) ~ (J -Z a )£r ]<p 2 + <p 3 

+[(nZ)*+(/-2*)f r ]^ - 4[/ + ±Z 2 J<P', ' 

in the region of the blades. The are correspondingly more 
simple for radial blades. However , the complexity of Q is pri- 
marily due to rotationaiity and compressibility. For example, if 
the fluid is incompressible and the inlet conditions are uniform the 
right side of Equation (64) becomes simply 

Q; = 2 U)r 2 j (114) 

regardless of the blade shape prescribed. 

-pecialized forms of the difference equations will appear in 
the various examples, o arts VIII. IX. and X. The general form of 
the difference operator is written as the sum: 

D 'M- f °;j <hj (115) 

the subscript i denotes the point of the net for which the operation 

is performed and the subscript j denotes the values of the stream- 
i unction which enter into the difference equation. Thus for the 
simplest difference operator, from Equation (110): 









I 
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D, (p) = <Pi, * (!-f r . j <P,-z i- <P,S +0*£) Pu - 4 Pic 

on = 0/3 = / ; a io = ~ 4- ( 116 ) 

O/z ~ l~TFi i a i* ~ '+TF,- 

The a ( j arc called "influence coefficients" in that they indicate the 

change in S s ^y , corresponding to a unit change in <p,-j ; or ~~ 

represents the change in tp/j due to a unit change in the "forcing 

function" &*/?■ whan all other are held constant. The a,, from 
fc/ tj J 

Equation (116) are shown for some point i: 




The difference equations must be satisfied at every point 
in the net. Tins is accomplished by outiuvell's relaxation tech- 
nique, whereby the streamfunction is modified until the desired 

♦ 

residuals are obtained. 

The solution of a finite difference problem is generally 
not the same as the solution of the corresponding differential equa- 
tion problem. I he amount of the discrepancy depends on the nature 
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of the function involved and on the sise of trie net used. There- 
fore it i c only necessary to obtain an accuracy in the solution of 
the difference problem commensurate with the discrepancy already 
involved . 

F) The Iteration Process 



rules which will govern the iteration process in all conceivable 
problems. It is, however, appropriate to discuss generally the 
two separate effects--compresaibility and vorticity. For the initial 
step of the iteration it is necessary to estimate the right sides of 
the equations. Consider, for example , the right side of Equation 



The manner in which rv( c/< ) may be estimated has already 
been discussed, in the examples solved the first estimates of 
rv( < p), so obtained, were accurate to within three percent of the 
final values (Figures 21 and 23) . The complete effect of rotation- 
ality, represented by the term in the bracket above, can likewise 
be expressed as a function of </> and r with comparable accuracy. 

In the relaxation process it is possible to improve the desired 
residual each time the streamfunction is improved, so that the 
relaxation and the iteration are actually performed simultaneously. 



it is not practical to attempt to formulate a general set of 



( 63 ) : 
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-eriodically rv(^ ) must be recalculated and a new, more accurate 
function of and r obtained. 

The compressibility effect cannot be handled so easily. If 
the Mach number is low and the rotationality effect is strong, as in 
the examples of Part IX, the first term, <fi n (/np) n , is much smaller 
than the last term and can be disregarded for the first estimation. 



In this case the effect of compressibility is not completely neglected, 
but appears in the f~J 2 i actor multiplying the rotationality term. 

It is true that under these circumstances the main effect of compres- 
sibility is it3 "influence" on the rotationality. For large subsonic 
h'ach numbers this approach is not possible. The first term may 
be of the same order as the second. In fact, in the examples of Part 
^ ’ where a constant tangential velocity was prescribed at the down- 



stream boundary, the two terms were exactly equal at the point on 
the downstream boundary where the axial velocity became sonic. 



i .da can be seen by writing equation (65), with oD = 0 and ij) ZI = 0, 



as 

r 7?(-F'Pr) - ' r r (Pw] ~ R 



it is well known that if the energy is uniform the local mass flow 
is a maximum when \v is sonic, hence/^ is aero when w i 3 sonic, 
mother viewpoint is to regard the residuals as forcing functions 
which act 10 displace the streamlines. 



i->ince the axial velocity can 



only be sonic at a ' throat", the forces must, act toward the sonic 
line so as to contract the streamlines there. j.Lc residuals, having 
a different sense on each side of the sonic line and being ’ continu-- 
ous", must therefore be sero on the sonic line. (hen v is not con- 
stant or when the energy is non-uniform, P w will not necessarily 
be maximum when w is sonic . The exact circumstances under 
which the residuals change sign ware derived in Part VI, where it 
was shewn that the residuals are sero on the line separating the 
"hyperbolic" and the "elliptic 1 regions, or, in other words, where 
the ‘governing velocity" is sonic. 

It is extremely important to use every means 'available to 
make the first estimate of the compressibility effect as good as 
possible. The rapidity of convergence depetids on the judgement 
used in the first estimation and each successive approximation. No 
general rules can be stated. In regions where the meridional velocity 
is nearly sonic, the streamsurfaces obtained during the iteration may 
be such that the channel between adjacent surfaces is "choked". 

When tnis happens the mass flow Pq s falls to the rignt of the density- 
curve in Figure 42 and the density is imaginary. If it is not possible 
to adjust the streamlines so as to have real densities, then it may be 
that ihe governing velocity is supersonic and the ' elliptic" form of 
the equation is not applicable, or even that no isentropic solution 
exists. If there is reason to believe isentropic solutions exist. 
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the '‘complete non-linear" equations derived in the appendix may 
be used for the transonic problem. Convergence is always slow 
when the governing velocity is nearly sonic, for the p vs pc^curve 
have nearly vertical tangents then. 
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This analysis of the flow in a mixed-flow compressor with 



blade shape prescribed is presented as an example of the three- 
dimensional motion of an incompressible fluid under the action of 
a system of rotating blades. Tirst a solution is obtained for the 
irrotational flow through the channel with no blades present. I'his 
serves as a basis of comparison for determining the additional 
velocities induced by the blades, and i.a also useful in making the 
first approximation to the rotational flow when blades are present. 



curved part of the channel boundary can be expressed analytically 
as 



where is the initial radius, r^ is the final radius, and L is the 
length of the curved cortion. The coordinates of the channel boun- 
do ries are given in Table I. 

In this example the total energy is uniform and the vorticity 
is zero at the inlet, station 0. Wore precisely, the tangential ve- 
locity is zero, the meridional velocity is axial, and the pressure 
is constant on the upstream boundary. In the blades and downstream 



I’he particular channel chosen is shown in figure 8. The 



z 



r- r , _ _£ 
r 2 - n ~ l 



s/n 2 n l 
2tt 



( 117 ) 



of the blades the total energy is non-uniform and the flow is rota- 
tional. The downstream boundary conditions are essentially these 
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t! at might Le expected if the ch nnel extended downstream to in- 
finity: the flov i3 independent of the axial coordinate . Thus the 
radial velocity is zero and the pressure gradient forces and the 
centrifugal forces are in balance on the downstream boundary. 

A) Irrotational Flow with no Blades Present 

For irrotational flow with no blades present the governing 
differential equation for the streamfunction (Equation (63)) is simply: 

*P tr ~ <Pr + 4* *2 z = O 0*8) 

and from Equation (UO) the corresponding difference equation is: 

* ('-§?)h + 4*3 i-i'+fA'P* - 4 <p a = O (119) 

The boundary conditions are: ip = constant on the hub and shroud 

contours t -yP r = tv- / on the upstream boundary, and -p *p z = u = o 
on the downstream boundary. Upstream p r is taken along the boun- 
dary, so that ip can be determined by simple integration and pre- 
scribed directly. Downstream the normal derivative is prescribed 
and p cannot be determined by integration along the boundary. The 
boundary value problem is therefore of the mixed type. However, 
since the flow far downstream is entirely independent of the axial 
coordinate all axial derivatives vanish and the differential equation 
(Equation 118) can be integrated easily. By this means p itself can 
be prescribed over all of the boundary. The mathematical problem 
for the irrotational flow is therefore relatively simple and consists 
of satisfying simultaneous!/ the difference equations for all points 
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in the net . 

fhe flow region is divided into a r.et or lattice and influence 
coefficients are calculated for each point. The net, the influence 
coefficients, and the boundary conditions are shown in .figure 9. 

The influence coefficients near the boundaries were calculated in 
the same way as Equations (106) and (108), but with a different 5 
for the "short legs" of the net. 

For a first approximation to ft the streamsurfaces were 
assumed to have the same shape as the boundaries. This is equiva- 
lent to assuming that the axial velocity is constant on each radial 
line. The actual residuals iff a for the smooth values of ft were 
calculated by means of the difference operator with a high degree 
of accuracy. These served as basic values of ft and (R a . The 
change in ft to make the residuals as small as possible was then 
determined by the relaxation process . As a final check residuals 
were again calculated by the difference operator. 

The final values cf the 3treamfunction and the corresponding 
actual residuals are shown in Figure 10. Upstream of station 12 
the smallest residual may be as large as 2 and downstream as large 
as 5 . This error corresponds to an error in the solution of the dif- 
ference problem of about one part in 2000. 

The axial and radial velocities for several representative 
stations are shown in Figures 11 and 12, respectively, the abscissa 



- 
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of Figure ll is the deviation from the mean axial velocity: 

l V - Wm 



& W = 






where w is the inlet velocity and w is the mean axial velocity. 

o ™ 



The ordinate is the non-dimen 3 icnal radius: 



where is the radius of the inner boundary and r Q is the radius 
of the outer boundary. 

The abscissa of Figure 12 is simply: 

~ u 

u = TvT 

I ines of equal pressure are shown in Figure 13 for flow 
through the channel with no blades present. The decrease in pres- 
sure occurring far downstream is due to the decrease of the cross- 
sectional area of the channel and the corresponding increase of the 
mean velocity. This should be taken into account when evaluating 
the pressure increase when the rotating blades are present. 

B) Rotatio nal Flow with Blades Present 

Rotational flow through the channel of the previous example 
is considered here . The vorticity is generated by a rotating sys- 
tem of blades which acts on the fluid between stations ( and 14. Ra 
dial blades are prescribed, hence the blade shape function f(r ,z) 
depends only on the axial coordinate z. based on the irrotational 
meridional velocities of the previous example the function f(s) was 
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chosen 30 that the rate of energy input, i.e. the rate of Change of 
the moment of angular momentum rv, was very small near the 
leading and trailing edges. The discontinuity in v across the lead- 
ing edges cannot generally be avoided, however the blades were 
selected such that the jump in v is positive. The rate of energy 
input is indicated by the slopes of the curves in Figure 20, where 
the desired and the actual ^ on the radius r = 2.5 are compared. 
The blade shape function is chosen to give the desired ~ on the 
basis of the irrotationai meridional velocities. The actual — is 
determined from the final meridional velocities with the prescribed 
blades present. The shape of the blades is shown in Figure 8 and 
the blade surface function is given in Table I. 

The applicable differential equations for the three regions 
of flow were derived in Part V-fi . 



whore rv = rv( <//) only and is evaluated at the trailing edge by Equa- 
tion (5C) , which here becomes: 



V Kr ~ J <Pr + 4*zz ~ & 



upstream of blades (66) 



(l + Z*) tp rr ~ (I-Z 2 lyip r + 



= 2 u)r Z in region of blades (o'* ) 



ip rr - -Lif/ r + p z ^ = (n/-U)r l ) downstream of blades (68) 




( 120 ) 
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The corresponding difference equations in regions where the net 
interval is § were given by Equations (110), (111), (113) and (114), 
which for the three regions become, respectively: 

h + ('-fMz + 4 a + to- o (i2i) 

P, + [(-§}]+ ( l+ b\ Z *] 'J'z + ^3 + [( , +2?\+( l ~fr)Z :! ] tp* 

-4[l + ^Z 2 ]<p o = 2cdrZS z ( 122 ) 

4>, * (l-Tr) 4>Z +-^3 + A * ft I h - 4 A ~ l™- Art 2 ( 1 23) 

where rv(<//) is evaluated by Equation (120). 

The influence coefficients for chese equations and fcr the 
equations which hold near the boundary where S is not constant are 
shown in Figure 14. 

The boundary value problem is the same as in Part VIII-.,. , 
except that here the streamfuncticn can net be prescribed exactly 
on the downstream boundary. Instead new downstream boundary 
conditions must be found for each approximation to rv(f). This is 
easily done by one-dimensional relaxation on the downstream boun- 
dary. The first approximation to the downstream boundary values 
was based on rv determined from the irrotational flow. For com- 
parison, the initial and final downstream boundary values are: 



r 


3.25 


3.0 


2.75 


2.5 


2.25 


2.0 


Initial </> 


0 


973 


1S5 7 


2654 


3367 


4000 


Final </, 


0 


977 


1863 


2660 


3371 


4000 
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iho p^ra.'iTetoro c. the problem, based on the maxin urn 
radius of 3 . 25 , are: 



Pressure coefficient: 


u)r * 
y 2 (u)r) z 


= 3.55 


Flow coefficient: 


W 

u>r 


=0.769 


Root-tip ratio: 


/© 
r i 


= 3 upstream 






= 1.625 downstream 




r 0 -n 
r i 


= 2 upstream 



= .625 downstream 



Blade aspect ratio: 0.232 

Static oressure increase: (-^ ) — (—] = 1.570 at tip, 

^fo / b/adms ' no biadoi 

- . 736 at root 

Le extreme tip-root ratio was chosen so as to amplify the 
three-dimensional properties of the flow, and blades of very low 
aspect ratio were orescribed in order to stress the importance of 
blade geometry. 

i’here was nothing unusual in the analysis. Relaxation and 
iteration were performed simultaneously, and the function rv( <p ) 
only had to be determined once after the first estimation, the 
initial and final values of r v(ip ) are plotted in Figure 21. The fin'l 
values of the strcamfunction and the corresponding desired and ac- 
tual residuals are shown in Figure 15 . The axial and radial veloci- 
ties are plotted fer several stations in .'igurcs 16 and 17 , using the 
same variables as already described. I ines of equal pressure are 
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si.own in figure i°. and lines of equal total energy in Tigure 19. 

The pressure was calculated by a relation derived from liquation 
(31): 

P~Po _ 2 u)rv f, + 7 

^PWo 2 ' I J 0 ^ 4 ) 



A comparison of the two solutions, with and without blades, 
does not lead to any extraordinary or unexpected revelations. 
There are, however, some interesting points that should be dis- 
cussed: 



i . The general effect of the blades on the meridional flo.v 
pattern is to increase the velocities near the shroud and 
decrease those near the hub. this may be considered as 
an improvement of the meridional flow, for the most likely 
place for separation is on the hub near station 13. Decreased 
velocities cm the hub therefore lessen the likelihood of sepa- 
ration . 

2. Ihe direction of the streamlines, given by , is changed 
very little by the action of the blades. 

3. The lowest static pressure, 7 ^ at station 13 on the 

s' Pm, a 

hub, is increased from -1.000 to -0.166 by the action of the 



blades . 



4. .although relatively large changes of axial velocity occur 
within the blade region, the downstream equilibrium values 




I 
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differ very little (4.6 percent) from the axial velocity 
with no blades present. However, the tangential velocity 
far downstream is approximately proportional to the ra- 
dius , and is of the same orcer as the axial velocity: 

~ - 1.146 at tip, - 7 - = 0.720 at root. Considerable 

diffusion would therefore be necessary to recover, as 
static pressure , all of the kinetic energy added. 

5 . Considerably more energy was added near the shroud 
than near the hub. This is unavoidable for radial blades, 
or, in fact, for any ' practical" blades of low aspect ratio. 

6 . The total velocity increases disccntinuously across 
the leading edge and it might seem that thi pressure 
should decrease correspondingly. There is, however, a 
sufficient increase of total energy to cause the pressure 
to increase. This is because the term co rv is always 
greater than the term -Lv 2 in Equation (124). 

7. The fact that the residuals are nearly constant on 
cylindrical surfaces extending downstream from the 
trailing edge of the blades meana essentially that the 
tangential vorticity is nearly constant on those surfaces. 
This confirms, for this particular example at least, the 
linearizing assumption used by h arble (rteferences 1 and 
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<i), when he assumed in the axial flow problem that the 
vcrticity was transported unchanged on cylindrical sur- 
faces . It is surprising that this is also true in this 
region , bounded by cylindrical surfaces but located 
downstream of a complicated region of mixed flow . 






31 - 



IX. AN c.K< L x'I : sV SUBSONIC ^LUv 
THROUGH N /CTUATOK C*I°-K 

.-.n example of the motion of an incompressible fluid acted 
upon by an infinite system of blades was presented in Part VIII. x's 
regards rotationality effects, compressible fluid motion through a 
system of blades can be analyzed similarly, the main difference 
being the density ratio multiplying the vorticity term. Furthermore, 
it was shown in Part V-_ that the flow in the region of the blades, 
where the vorticity is bound, is much less complicated than in the 
region downstream of the blades, when the vorticity is free. There- 
fore the present example, conceived in order to isolate the role 
compressibility plays in altering the rotationality effects, is con- 
cerned only with the more complicated upstream and downstream 
regions . f he blade region is therefore concentrated into an "ac- 
tuator disk" in which the tangential velocity jumps discontinuousiy 
from zero to some finite value. It is assumed that the "blade sys- 
tem" is stationary, hence the total energy of the fluid is constant 
throughout the field. The matching conditions for the meridional 
flow at the discontinuity are the same as stated in Part V-D for 
leading edge discontinuities . 

The channel boundaries are concentric cylindrical surfaces , 
and at the entrance the axial velocity is constant and the radial and 
tangential velocities are zero. The jump in tangential velocity is 
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e^ual in magnitude to the inlet velocity of the incompressible example, 
.'or the downstream boundary condition the flow is independent of the 
t-xial coordinate, ibis problem is solved for two values of tbe inlet 
teach number (0 and 0.2), the inlet mass- velocity and the tangential 
velocity being the same in both cases. 

The applicable equations, obtained from Equations (63) ami 
(65) with tO = 0, (rv)| =0, and = constant, are: 

4>rr - -p (p r + <piz - (lnP) n upstream 



= <p n (foF) n -±(?)* £L(£ji) z downstream 



( 125 ) 



where rv(^) and — (p~- ) z are functions of d> only and are evaluated 

dp ' 7 

at the downstream side of the discontinuity. 

~he boundary conditions are: ip - constant on the cylindri- 

cal boundaries, -~p r = 1 on the upstream boundary, and ^ = 0 
on the downstream boundary. For each approximation to rv(<p) 
the downstream boundary condition can be modified so that p is pre- 
scribed, as in the previous example (Part VIII-B). 

A) Incompressibl e Flow Through the Actuator 

The difference equation corresponding to Equation 125 for 
incompressible flow is: 

^('-£1^2 <-^3 -4^0 = O upstream 



(126) 



Lsi{rv ) z c 2 downstream 

~ 2 dp ° 



The inlet velocity is taken as 1000 and the boundary conditions are : 
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- 500 (9 - r ) on upstream boundary 
p ~ 0 on r = 3 
tp - 4000 on r = 1 
^ = 0 on downstream boundary 

The influence coefficients and the net points of the flow field are 
shown to scale in Figure 22. 

The first approximation to rv( p ) was based on irrotational 
flow within the same boundary. That is, it was assumed that the 
upstream value of p held throughout the field. The rv so obtained 
is shown in Figure 23. Before relaxation, the corresponding down- 
stream boundary values were further improved by assuming half 
the change in <p occurred upstream of the actuator , thus obtaining 
an approximate value of <p at the actuator and consequently a new 
rv(</> ). The boundary values corresponding to this last rv were 
then used for the first complete relaxation, after which rv( p) was 
found to be unchanged. The result of this complete relaxation was 
therefore the final solution of the problem. 

Values of the axial, radial, and tangential velocities at sev- 
eral representative stations are shown in Figures 24 , 25, and 26. 
Lines of constant static pressure are shown in Figure 27. Far down- 
stream of the actuator the axial velocity is much greater near the 
hub than near the shroud. 7’he maximum value of the radial velocity 
occurs ju3 1 beiow the "center' of the channel. The pressure is 
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discontinuous acres a the actuator, and most of the pressure changes 
occur downstream of the actuator . 

Perhaps the most significant feature of the solution is that 
the radial velocity and the changes in axial velocity are approximately 
symmetrical about the actuator. It is clear, however, that they arc 
not exactly symmetrical since the forcing function is not symmetri- 
cal. One consequence of linearisation (Reference 4) is that the mer- 
idional flow is symmetrical about the discontinuity. 

The solution cf this example will be discussed and compared 
with the compressible solution of the next section. 

B) Compressible Flow Through the /iCtuator-V-j n j e ^ = 0.2 

The difference equation corresponding to Equation (125) is: 

P, - 0-fr)^ + <P 3 +("-£f)^ - + 

upstream (127)* 



" frn Mn ~ i ] & 2 downstream 

The boundary conditions for the streamfunction* were essentially 
the same as for the previous example (IX~->). On the upstream 
boundary <P r ~ w - 1000, and - 0.2. Cther parameters on the 



upstream boundary were: 



p 

Po ~ 




- 1 .020 


a 

do 


(l-tirM*)* 


= 1 .004 


W - 


/ooo & 


= 1020 



’ .. or this exa-.mle only, the streamfunction is defined by: 
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=5100 

a 0 = 5080 

The tangential velocity jumped diacontinuou.sly to 1000 in the actuator. 

The influence coefficients and the net points of the field are 
the same as in the previous example (Figure 22). 

The first step in obtaining a solution was to estimate as ac- 
curately as possible, by any means whatever, the right side of the 
equation. This was done as follows: The vorticity term, 
was very closely approximated by simply using the value from the 
incompressible solution iust obtained. The extreme accuracy of 
this choice is indicated in Figure 28, where the initial and final val- 
ues of r v(ip ) are plotted. The first estimation of the density term 
was likewise based on the irrotational velocities, but only values 
on the downstream boundary were calculated first. When this was 
done it was clear that the streamlines deflected les3 when compres- 
sibility was present, tor the compressibility term, (lr>f} n 5 1 » was 
negative and subtracted from the vorticity term, £- U- (rr) z J 

which was positive. In addition (~f~) Z > which decreased according 
to the increase in total velocity downstream of the actuator, also 
acted to reduce the right side of the equation. However, the right 
side of the equation was estimated on the downstream boundary using 
the density term just calculated, and the streamfunction and the cor- 
responding density were calculated there. These values were treated 
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cio downstream boundary values for the first complete relaxation, 
in order to complete the estimation of the residuals throughout the 
field the density changes calculated from the irrotational stream- 
lines were corrected by a factor chosen to give the correct density 
on the downstream boundary. In this way it was possible to make 
a rather accurate estimate of the right side of the equation, and 
the desired residuals for the first complete relaxation, From the 
results of the first relaxation the densities and velocities, and 
finally new desired residuals, were obtained. The new desired 
residuals agreed very well with the previous ones except in regions 
where the tv.ach number was greatest. Small changes of the stream- 
function brought the residuals into satisfactory agreement and led 
to the final solution, presented in Figures 23 through 32. 

i'he velocities are expressed in terms of the inlet axial 
velocity, rather than a velocity of sound, for comparison with the 
incompressible solution. The h-'ach number based on the total ve- 
locity is shown in Figure 32. I’he velocity components are similar 
to those of the incompressible solution, except of course for the 
discontinuities at the actuator. 

C) .. Discussion of the Solutions 

In the two foregoing examples an otherwise axial flow is 
distorted by vorticity generated by a so-called actuator disk. The 
strength of the vorticity and the boundary conditions of the flow are 
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tuc s>umo in beta examples , yet the solution s , although generally 
similar, differ distinctly in at least two way;,: the total deflection 
of the flow brought about by the vorticity is decidedly less when 
the flow is compressible; the axial variations of the compressible 
flow are concentrated more in the vicinity of the actuator, the 
concentration being more intense where the ^ sch number is higher, 
’’he first point of difference was discussed in the previous 
section and is confirmed by Figure 33, where the streamlines of 
the two solutions are compared. A general conclusion set forth 
in Part VI is that the locai effect of vorticity, as regards defiec- 
ticn of tne strearosurfaces , reduces to zero as the "governing" 

Vach number increases to one, and then increases, but with oppo- 
site sense, as the Vlach number increases above one. The second 
point of difference could probably be predicted by the Prandtl- 
t/lauort similarity transformation, applied to linearized equations. 
Its occurrence here is evident in Figure 32, and to some extent 
in r igures 25 and 30 where the maximum of the radial volocity 
for the compressible flow occurs nearer the center of the channel, 
indicating more rapid axial changes near the hub where the 5>/ach 
number is larger. 
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X. t.-'l-L CF TR-vNSGNIC JLC.. 

r K#i r ugi: an ac ru. > for di-i: 

This section provides examples of flow through actuator 
disks with maximum JS/ach numbers near one. Two examples , 
with the same inlet Vach number (0.555), but with actuator disks 
Oi different strength, are considered. In one example the maximum 
h/ach number attained is one, in the other, 1.12. The actuator 
model is the same as described in Part IX but the tangential ve- 
locity induced at the actuator is not constant there . Instead the 
tangential velocity is prescribed on the downstream boundary, and 
the jump at the actuator is not known until the problem is solved. 
Aith this prescription it is possible to replace the downstream 
boundary condition, || =0, with the condition that the stream- 
function itself is prescribed, for the downstream conditions can 
be determined without knowing the intermediate flow. The limit- 
ing flow far downstream of the blades was discussed in Part V-F . 
j or tue particular case wnere the tangential velocity is constant, 
Figures 43 and 44 show the maximum Ivfach number as a function 
of the inlet l,,ach number and the tangential velocity. It was also 
shown in Part V-F that there is a maximum tangential velocity’ 
ti.at can be imposed by an actuator disk and that this maximum 
corresponds to the choking condition. The limiting conditions 
for these examples are indicated by the two points shown in Figure 



43 . 
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upstream 

downstream 



( 128 ) 



ihc differential equation (Equation 125)) of the previous 
examples is applicable here but the definition of the streamfunc- 

tion * is ali 8 htl V different. The difference equation for both ex- 
amples is (from Equation (111)): 

+ ('-£ \K + ^3 + (t+ <l>4~ 4-lp 0 

“ [^n (In f)n ] £ 2 

-MM, -rfrS^yjs * 

•) Flow with Maximum l/lach Number of C ne 

The following conditions are prescribed: 

Inlet E.ach number: lv! = .553 , b/l'* = 58"^ 

Tangential velocity: v = 0 at inlet 

(~a$~ *21 on downstream boundary 

and calculated from these are: 

Maximum axial velocity: — = .888 

Maximum ivfach number: 1 

The influence coefficients for each point in the net and the boundary 
values for p are shown in Figure 22. 

he method of solution was essentially the same as de- 
scribed in the previous example. In determining the approximate 
residuals for the first relaxation the value of the rotationality term 
2^® obtained on the downstream boundary was assumed constant 
* r or tlxese examples the streamfunction p is defined as: 



— = -&--L d) , • iv _ jf 

P ^ P r ^ r 



a+ 



P r 
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on cylindrical surfaces extending upstream to the actuator . The 
density throughout the field was estimated using the known values 
of the density upstream and downstream and the approximate mag- 
nitude of the density discontinuity at the actuator. The first de- 
sired residuals were obtained in this way and a complete relaxa- 
tion performed. The final solution then followed simple iteration 
in which the desired residuals were obtained from the solution of 
the previous step, lifter tv o complete relaxations the residuals 
were in good agreement except in the region where the h/ach num- 
ber was largest. >/inor modification then led to the final solution. 

Che axial and radial velocities are shown for representa- 
tive stations in Figures 35 and 36. The Iv.ach number based on 
the total velocity is shown in Figure 37. 

This solution exhibits the properties already discussed in 
the previous compressible example of °art IX, as noted in the 
comparison o<. the next section. 

B) Flow with h axitnum h'ach Number of 1 . Id 

vhe following conditions are prescribed: 

Inlet 2- ach number: M = .555, V** = .584 
Tangential velocity v = 0 at inlet 

(~^~J - -264 on downstream boundary 

and calculated from these are: 
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Vaximum axial velocity = - .966 

Maximum axial Vach number = — = .935 

a 

» S ~t 

Maximum total Vach number = [ w J = 1.12 
The influence coefficients for each point of the net and the boundary 
values for ip are shown in Figure 38. In addition the velocity, den- 
sity, etc. , at the downstream boundary, are recorded in Table II. 

The residuals were first estimated using the meridional 
streamlines of the previous solution. A complete relaxation was 
performed and the iteration started. Unfortunately the densities 
obtained from the second complete relaxation were imaginary and 
straightforward iteration could not be continued. To remedy this 
the axial mass- velocity , , was changed just enough to make 

the densities real. Then the iteration could be continued. It was 
observed that the numerical solution of this example was less ac- 
curate than the solution of the previous example , although the net 
points (Figures 34 and 38) were four times as dense. This is due 
to the nearly' vertical tangent of the density curve (Figure 42) when 
the m;ridional velocity' is nearly sonic, and to the fact that the ra- 
dial derivative of the streamfunction could not be accurately deter- 
mined on the inner radius where the large Vach numbers occur. 

The axial and radial velocities are shown in Figures 39 and 
40, the Vach numbers in Figure 41. The solution exhibits the same 
properties as have already been described in the preceding examples. 
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It is surprising that the streamsurfaces are deflected less in this 
second example where the tangential velocity is greater. This can 
be seen by comparing the downs tr earn values ol the streamiunction 
in figures 34 and 33, and indicates that there is a certain tangen- 
tial deflection for which the deflection of the streamsurfaces is 



maximum . 




I 






I 






- u- 



*:i. cl nci iniotr 

theory has been developed which permits the analysis of 
compressible flow in turbomachines having infinitely many blades 
when the governing velocity is subsonic. Several examples, solved 
by the method of finite differences, have been presented. The funda 
mental idea underlying the theory is that the force field represent- 
ing the infinitely many blades is necessarily a "pseudo-conservative 
fieid. Because of this the three components of the field can be ex- 
pressed in terms of two functions, one describing the input of 
energy, the other the shape of the blades. The functions which 
must be prescribed and the boundary conditions which must be im- 
posed are then quite clear, and the heretofore more difficult direct 
problem becomes relatively easy. 

Two principal effects of compressibility were noted: the 
deflection of the streamsurfaces brought about by a given vorticity 
distribution decreased as the Mach numbers (subsonic) were in- 
creased; the streamwise variations of the flow became more con- 
centrated as the ..ach numbers were increased. 

[he foregoing theory is limited in that the fluid must be 
non-viscous and the number of blades must be infinite. In addition 
the application of the theory when the "governing" velocity is super- 
sonic is questionable. It is clear that non-viscous flow under the 
action of a finite number of blades must be completely understood 
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before the vUcous, problem is. attacked. Viscosity has no place in 
problems in which the number of bladej is infinite, for the blades 
do not act as boundaries of the flow. 

ihe e:iton " ion io m< chines with a finite number of blade* 
could conceivably be carried out by three-dimensional finite differ- 
ence methods but the analysis would be very lengthy. A first order 
approximation to the pressure and total energy is indicated in the 
appendix. Perhaps the least understood problems are those in which 
the governing velocity is transonic. The equations were developed 
and discussed generally in Part VI. /. suggested approach might 
be to consider onLy the region in and upstream of the blades, where 

thC V ° rticity effects are linear * Prescription of radial blades 
would allow further simplification by means of the transformation 
given in liquation (121). Then a perturbation equation analogous to 
the transonic equations could be developed. Any analysis which 

" lH lead to a better understanding of the cushioning effect men- 
tioned in Fart VI would be a worthy contribution. 
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- .i f»SNDIX 

jj; ir a t ( rder Approxi mation to 

the lo v with a Finite Number of blades 



F first order approximation to the flow with a finite number 
of blades can be obtained from the solution with an infinite number 
of blades by replacing the body force field acting in a sector between 
two blades with a pressure gradient force, thus essentially revers- 
ing the reasoning which first led to the assumption of infinite blades. 



Then integration of the pressure gradient from one blade to another 
will give the circumferential pressure distribution between the blades. 

The mean value of the pressure, as obtained from the anal- 
ysis with an infinite number of blades, is denoted by p . Let p' 

denote the deviation of the total pressure from the mean: p = p -‘■■o' 

1 m 1 

* he force field is then replaced by the pressure gradients as follows: 
Fr ~ ~ A -f r - 



/ Up’ 
P ar 



rr _ i / _ _ / <>*>' 

i9 A r ~ ~ r a& 

fr — -J r - - JL 

Fz ~ - f* 02. 

If there are n blades distributed uniformly around the axis , the 
second equation can be integrated from one blade to the next to give 
the pressure discontinuity across a blade: 

because A is independent of ■& , the tangential variation of the pres- 
sure is linear between blades, /he deviation from the mean pressure 



can therefore *jc eusiiy determined t all noints in the region of the 
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blades. The corresponding approximation of the total energy, oO rv, 
can be obtained from a form of Equation (31), using the velocities of 
the solution with an infinite number of blades. 

it is significant that the equilibrium of the fluid has not been 
disturbed by the introduction of the pressure force for the force fie la, 
nor has fulfillment of the continuity requirement been altered. Ac- 
tually the only assumption involved is that the tangential variation 
of the pressure is linear, i.e., X i3 independent of O . This first 
order approximation therefore seems very reasonable, particularly 
•when there are many blades, and may be satisfactory in other cases. 
Unf o r tuna te 1 y exact solutions with a finite number of blades are not 
available for determining the accuracy of the approximation. 
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HUB AND SHROUD CONTOURS- BLADE SHAPE 
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FIG. II 

AXIAL VELOCITY- CLEAR CHANNEL 
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RADIAL VELOCITY- CLEAR CHANNEL 
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The Channel Boundaries and the Blade Shape Parameters for 
the Mixed Flow Example 
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Limiting Values Far Downstream for Transonic 
Actuator Example with M max = * • * ^ 
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Tabulated Values of the Isentropic Energy delation 
Density ?.atio vs. I asa Flow and effective Tangential Velocity 
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